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Abstract. Medium resolution (A/A.A = 5 000 - 10000) VLT-ISAAC M-band spectra are presented of 39 young stellar objects 
in nearby low-mass star forming clouds showing the 4.67 pm stretching vibration mode of solid CO. By taking advantage of the 
unprecedentedly large sample, high S/N ratio and high spectral resolution, similarities in the ice profiles from source to source 
are identified. It is found that excellent fits to all the spectra can be obtained using a phenomenological decomposition of the 
CO stretching vibration profile at 4.67 pm into 3 components, centered on 2143.7 cm"', 2139.9 cirT 1 and 2136.5 cm"' with 
fixed widths of 3.0, 3.5 and 10.6 cm" 1 , respectively. All observed interstellar CO profiles can thus be uniquely described by a 
model depending on only 3 linear fit parameters, indicating that a maximum of 3 specific molecular environments of solid CO 
exist under astrophysical conditions. A simple physical model of the CO ice is presented, which shows that the 2139.9 cm" 1 
component is indistinguishable from pure CO ice. It is concluded, that in the majority of the observed lines of sight, 60-90% of 
the CO is in a nearly pure form. In the same model the 2143.7 cm" 1 component can possibly be explained by the longitudinal 
optical (LO) component of the vibrational transition in pure crystalline CO ice which appears when the background source 
is linearly polarised. The model therefore predicts the polarisation fraction at 4.67 pm, which can be confirmed by imaging 
polarimetry. The 2152 cm" 1 feature characteristic of CO on or in an unprocessed water matrix is not detected toward any source 
and stringent upper limits are given. When this is taken into account, the 2136.5 cm" 1 component is not consistent with the 
available water-rich laboratory mixtures and we suggest that the carrier is not yet fully understood. A shallow absorption band 
centered between 2165 cm" 1 and 2180 cm" 1 is detected towards 30 sources. For low-mass stars, this band is correlated with 
the CO component at 2136.5 cm" 1 , suggesting the presence of a carrier different from XCN at 2175 cm" 1 . Furthermore the 
absorption band from solid 13 CO at 2092 cm" 1 is detected towards IRS 51 in the p Ophiuchi cloud complex and an isotopic 
ratio of 12 CO/ 13 CO = 68 ± 10 is derived. It is shown that all the observed solid 12 CO profiles, along with the solid 13 CO profile, 
are consistent with grains with an irregularly shaped CO ice mantle simulated by a Continuous Distribution of Ellipsoids (CDE), 
but inconsistent with the commonly used models of spherical grains in the Rayleigh limit. 
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1. Introduction 

Solid interstellar CO was first reported bv ISoifer etaD Jl979h 
who observed a strong absorption band at 4.61 fim in a spec- 
trum from the Kuiper Airborne Observatory (KAO) along the 
line of sight towards the massive young stellar object (YSO) W 
33A. Although the absorption band later turned out to be domi- 
nated by a CN-stretching carrier, this was the first observational 
indication that the ice mantles on dust grains in dense clouds 
carr y a rich chemistry, as in dicated by early chemical models 
fe.g. lTielens & H agen 1982, and references therein). Until this 
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time the early ideas about "dirty ice" were mostly supported 
by observations of the 3.08 /im water band and the extended 
red wing of this band, which was sug gested to be due t o NH? 
- H?P complexes and CH 3 OH (e.g. iGillett & Forrestl fl973l 
IWillneretallll982l) . The presence of both solid C O and XCN 
toward s massive young stars was confirmed by lLacv et all 
The identification and subsequent study of the ices in 
general, and of the CO stretching mode at 4.67 /jm in particular, 
was aided by laboratory spectroscopy of ices inten ded to simu- 
late in t erstellar ices, such as th o se presented in e.g. | Hagen et alJ 
dl979h : ISandford et al.1 dl988l): ISchmitt et all dl989l) and m ore 
recently by lEhrenfreund et all d 19971) : lElsila et al.1 i 19971) and 
iBaratta & Palumba Jl998l) . The typical spectral resolution of 
these studies was A/AA ~ 1 000 - 2000 (1.0 - 2.0 cm" 1 ). 
However, at this resolving power the I2 CO stretching mode as 
well as the I3 CO stretching mode at 2092 cirT 1 are in general 
not fully resolved. This problem has caused difficulties when 
interpreting both interstellar and laboratory spectra. 

The laboratory studies have shown that some band profiles, 
and in particular that of the CO stretching vibration mode, are 
sensitive to the presence of secondary species in the ice as 
well as to thermal and energetic processing of the ices, and 
attempts have been made to quantify some of these effects. For 
instance, the 4.67 fim CO stretching vibration profile has been 
used in attempts to constrain the concentration of species like 

CO?, O?, No and the deg ree to which the CO is mixed with the 

■ ' 11 1 1 11 1 1 — 1 

water component (Tielens et al. 1991; Chiaret al. 1994, 1995, 

Il998h . This informatio n is otherwise hard or impos sible to ex- 
tract from observations ( Vandenb ussche et alJl999h . It has also 
been suggested that the CO profile may be used as an indicator 
of the degree to which the ice has been thermally processed if 
different components of the profile correspond to CO ices in en- 
vironments of different volatility. For instance, pure CO may be 
expect ed to desorb at lower temperatur e s than CO mixed wit h 
water dSandford & Allamandolalll988t ICollings et aljEool). 
The laboratory simulations have been used for comparisons 
with a large number of spectra from both high and low mass 
young stars and quiescent cloud lines of sight probed towards 
background late-type giant stars, to constrain the ice mantle 
composition for a wide range of different cloud conditions. 
Often these studies employ the full range of available labora- 
tory spectra in order to fi nd the best-fitting combination of ice 
mixtures to each source (|Tjejens_e^alj|lj9j| jKerret al.ll 19931 
IChiar etalJll994l Il995l Il998t iTeixeira et alJll998h .~This ap- 
proach is known to suffer from serious degeneracies and is fur- 
ther complicated by the fact that profiles from ice mixtures with 
a high concentratio n of CO are significantly affected by the 
shape of the grains (Tie lens et all 199 II) . Consequently, the de- 
tailed environment of solid CO in interstellar grain mantles and 
its role in the solid state chemistry is not strongly constrained 
by observations. Clearly, a large, consistent sample of lines of 
sight is necessary to obtain further information on the "typical" 
structure of ices in space. 

Recent results from medium to high resolution (R - 5 000- 
25 000) spectrometers mounted on 8 m class telescopes are 
beginning to cast new light on the 3-5 yum region of proto- 
stellar and interstellar spectra. In particular, new insights have 
been gained through careful studies of the resolved CO ice 



band s tructure and its relation to gas phase CO dBoogert et alJ 
20021)1) , using the hig her sensitivity to probe new environ- 
ments llThi et alJl2002l) and using the 13 CO ban d to constrain 
both c omposition and grain shape corrections dBoogert et alJ 
l2002al) . 

Here we present medium resolution (R = 5 000 - 10000, 
Av = 0.4 - 0.2 cirT 1 ) M-band spectra of 39 low mass young 
stellar objects obtained with the Infrared Spectrometer And 
Array Camera (ISAAC) mounted on the Very Large Telescope 
(VLT) ANTU telescope of the European Southern Observatory 
at Cerro Paranal. These spectra have allowed us to obtain mean- 
ingful statistics on the general shape of the CO ice profile for 
the first time. In this work, we show that statistically there is 
enormous advantage in using a large sample of fully resolved 
spectra from a well-defined set of young stellar objects. By us- 
ing the highest possible spectral resolution, we are able to cast 
new light on the underlying ice structure from a detailed anal- 
ysis of the CO stretching vibration profile. The L-band spectra 
and the analysis of the gas phase CO lines seen in the M-band 
spectra will be presented in later articles. 

This article is organised as follows. In Sec.[5]the sample 
of sources is described including the adopted selection criteria, 
observational procedures and main features in the spectra. Sec. 
|3] presents a phenomenological decomposition of the CO ice 
band and the observational results of the decompositions are 
then discussed. A simple physical model is described in Sec. 
0] which can partly explain the success of the phenomenolog- 
ical decomposition. The results from the observations and the 
modeling are discussed in Sec. |5] and the astrophysical impli- 
cations are outlined. Finally, Sec. |6] presents the conclusions 
and suggests further observations and laboratory experiments 
which can be used to test the results presented in this work. 
Appendix|A]reviews the simplest mechanisms governing solid 
state line shapes and AppendixlBlgives comments on individual 
sources in the sample. 

2. The sample 

2.1. Source selection 

The primary objective was to perform a broad survey of the 
dominant ice components on grains surrounding low-mass em- 
bedded objects in the nearest (D < 500 pc) star forming clouds 
in the southern sky. Emphasis was put on obtaining good S/N 
ratios at the highest possible spectral resolution in order to 
study not just the main H2O and CO features but also the 
fainter ice species, probe sub-structure in the ice features, and 
obtain a census of gas-phase ro-vibrational emission and ab- 
sorption towards low mass young stars. The observed sample 
was primarily extracted from the list of sources to be observed 
with the Space Infrared Telescope Facility (SIRTF) as part of 
the L egacy project From Molecular Cores to Planets (lEvansI 
l2003h . but was supplemented with additional sources. The ma- 
jority of sources have luminosities of Lboi < 10 L , with only a 
few brighter sources being included for comparison. Embedded 
sources are selected from their published Lada classes if no 
other information is available indicating youth. The sample in- 
cludes class I sources as well as flat spectrum sources defined 
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by having a vanishing spectral index between 2.5 and 15 fim 
(Bon temps et alJl2001l) . In addition, the background star CK 2 
in Serpens was observed to probe the quiescent cloud material. 

A number of the sources have published M-band spectra, 
albeit at lower spectral resolution, coverage and S/N ratio and 
were re-observed in order to get a consistent sample. The list 
of observed sources and relevant observational parameters is 
given in Tabled 

2.2. Observations and data reduction 

The ISAAC spectrometer is equipped with a 1024 x 1024 
Aladdin array, which is sensitive in the 1-5 fim region. The 
pixel scale is 0.148" pixel -1 when using the spectroscopic ob- 
jective. The telescope was operated using a standard chop-nod 
scheme with typical chop-throws of 10-20". Both L and M- 
band spectra of each source were obtained. The medium reso- 
lution grating was used for the 4.53 -4.90 fim region, resulting 
in resolving powers of A/AA = 5 000 - 10 000 (0.4 - 0.2 cirT 1 ). 
Here a single setting covering 4.53 to 4.77 /mi was always ob- 
served and a second setting extending to 4.85 - 4.90 fim was 
obtained if necessary, e.g. because of the presence of strong gas 
phase lines in the first setting. Often, it was possible to place a 
secondary source in the slit by rotating the camera. In this way 
spectra of both components of all binaries with separations of 
more than ~ 1" were obtained. 

The data were reduced using our own IDL routines. The 
individual frames were corrected for the non-linearity of the 
Aladdin detector, distortion corrected using a startrace map, 
and bad pixels and cosmic ray hits were removed before co- 
adding. Both positive and negative spectral traces were ex- 
tracted and added before division by the standard star spectrum. 
An optimal small shift and exponential airmass correction be- 
tween the source and the standard were applied by requiring 
that the pixel-to-pixel noise on the continuum of the final spec- 
trum is minimized. Before an exponential airmass correction 
is applied, the detector and filter response curve has to be re- 
moved in order to obtain a spectrum of the true atmospheric ab- 
sorption. The combined filter and detector response curve was 
obtained by fitting an envelope to both the standard star spec- 
trum and a spectrum of the approximate atmospheric transmis- 
sion and taking the ratio of the two envelopes. 

The spectra were flux calibrated relative to the standard and 
wavelength calibrated relative to telluric absorption lines in the 
standard star spectrum. The uncertainty in the flux calibration 
is estimated to be better than 30% for most sources and the 
wavelength calibration is accurate to ~ 5 km s . The high ac- 
curacy in wavelength calibration is due to the large number of 
well-defined and separated telluric lines present in a medium 
resolution M-band spectrum. For some of the weaker sources, 
the wavelength calibration becomes considerably more inaccu- 
rate because of difficulties in determining the small shift always 
present between the standard star spectrum and the source spec- 
trum. 

Finally, optical depth spectra were derived by fitting a 
blackbody continuum to the regions where no features are ex- 
pected to appear, namely around 4.52 - 4.55 fim and around 



4.76 - 4.80 fim. In every case care was taken not to include gas 
phase lines and regions affected by strong telluric residual in 
the fit. In cases of doubt the fit was compared to the continuum 
in the 4.0 - 4.2 yum region and adjusted accordingly. However, 
normally the L-band spectra were not used for the continuum 
fit, due to uncertainties in the relative flux calibrations between 
the two spectra. 

2.3. Features in the M band spectra 

The most prominent features in the 4.5 - 5.0 fim region of low 
mass young stellar objects are the absorptions from the stretch- 
ing mode of solid CO around 4.67 fim and the multitude of 
lines from the ro-vibrational fundamental band of gas-phase 
CO extending across the entire M-band. In addition the re- 
gion covers the CN stretching band around 4.61 /mi as well 
as the bright Pf/3 hydrogen recombination line at 4.65 fim and 
the S(9) line from molecular hydrogen at 4.687 fim. The VLT 
spectra in general show the presence of all these features, al- 
though the shape and strength of each feature varies greatly 
from source to source. For example, the ro-vibrational lines 
of CO are seen as both broad and narrow absorption lines 
and as broad and narrow emission lines. The broad absorption 
lines are often asymmetric indicating the presence of outflows 
and infall. The broad emission lines show a higher degree of 
symmetry and are clearly caused by hot gas as indicated by 
the presence of CO molecules excited to the second and even 
third vibrational levels. One source showin g narrow emission 
fines ( GSS 30 IRS 1) is treated in detail in Pontopp idan et alJ 
All the gas-phase CO lines will be analysed in a later 
paper (Pontoppidan et al., in prep). At a resolving power of 
A/AA = 10000 (30 km s" 1 = 0.2 cnr 1 ) all solid state features 
are completely resolved, including the central parts of the main 
CO feature and the very narrow 13 CO feature. This is generally 
not the case for previously obtained astronomical and labora- 
tory spectra at resolutions of A/AA = 1 000-2 000 (1 -2 cm -1 ). 

The individual M-band spectra with detected CO ice are 
shown in Figs.^Qon an optical depth scale. IRS 46, IRS 54, 
LLN 39 and LLN 47 have no or only marginally detected CO 
ice bands an d are not shown. The spec trum of GSS 30 IRS 1 is 
presented in Pontoppid an et alJ (12002) and shows the presence 
of a very shallow C O ice band. Th e spectrum of CRBR 2422.8- 
3423 is shown in iThi et alJ ll2002h and the s pectrum of LLN 
17=IRAS 08448-4343 is shown in lThil d2002h . 

2.4. New detection of solid 13 CO 

The 2092 cirT 1 band of solid I3 CO has been detected in one 
source in the sample, namely IRS 5 1 in the p Ophiuchi cloud. 
The source is a transitional Class I-II object and shows a very 
deep (t ~ 5) I2 CO ice band. The 13 CO band has an optical 
depth of 0. 12 and a width of 1 .2 crrT 1 . The detection is analysed 
in Sec. 14.41 This is only the secon d 13 CO feature repor ted in the 
literature after NGC 7538 IRS 9 (Boo gert et all2002al) and the 
first towards a low mass source. 
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Fig. 1. VLT-ISAAC M-band spectra of sources in the p Ophiuchi cloud. The spectra have been put on an optical depth scale 
by fitting a blackbody continuum to parts of the spectrum free of intrinsic features as well as deep telluric lines. In general, 
the continuum points are taken in the region 4.52^4.55 //m and between 4.76^.80 fim. The thick solid curve shows the total 
phenomenological model fit to the spectrum and the thin solid, dashed and dot-dashed curves show the individual components 
(blue, middle and red, respectively). The inset figures show the fit in the region of the 2175 cm -1 feature. 



3. Analytical fits 

3.1. A phenomenological approach 

One approach to the analysis of interstellar ice bands, es- 
pecially in the case of the CO stretching vibration, is via a 
mix-and-match scheme, whereby laboratory spectra of differ- 
ing chemical composition, relative concentrations, temperature 
and degree of processing are compar ed to astronomical spec- 
tra to find the bes t fits (e.g. IChiar et alJ Il99l 1 19951 Il998t 
iTeixeira et al.ll998l) . While the presence of solid CO can be un- 



ambiguously identified using the laboratory data, the mix-and- 
match approach is seriously affected by the degeneracy intro- 
duced by the dependancy of the CO vibration profile on many 
experimental parameters. In Sec. l5.3l we will return to this point 
and discuss the best strategies for comparing laboratory data to 
interstellar spectra. 

Due to the degeneracy encountered when attempting to 
match CO ice bands in space with a database of interstellar ice 
simulations, a more phenomenological approach is attempted 
here. Results from laboratory work can then be applied sub- 
sequently to explain the general trends identified using a phe- 
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Fig. 2. p Ophiuchi sources, cont. Legend as in Fig.[0 



nomenological method of analysis. Similar approac hes to the 
analysis of ice bands were used by Keane et ail il200lh fo r the 
6.85 ju m band. lBoogert et alJ J2002al) for the CO band andFrhl 
(2002) for the 3.08 /im water band, although for smaller sam- 
ples or for individual sources. 

The motivation is to drastically reduce the number of free 
model parameters and at the same time break the degeneracy 
by requiring the ice profiles of the entire sample to be de- 
scribed by a single model. The literature suggests that 3 dis- 
tinct components are fundamental to the CO stretching vibra- 
tion profile. It is well known that the CO band can be decom- 
posed into a narrow component centered around 2140 cm -1 
and a broader red-shifted component situated somewhere in the 
region 2130 — 2138 c m" 1 dTielens et al.l ll 991). Furthermore, 
Boo gert et al.1 d2002bl) detected the presence of a third narrow 
blue component at 2143 cirT 1 in the spectrum of the low mass 
young star L1489 in Taurus. We will in the remainder of the 



text refer to the 2134 cm feature as the red component (rc), 
the 2140 cirT 1 feature as the middle component (mc) and the 
2143 crrT 1 feature as the blue component (be). 

The qualitative identification proposed in the literature of 
the three components is as follows: The red component is tra- 
ditionally identified as CO ice in a mixture with hydrogenated 
molecules such as H2O or NH3, while the middle component 
is thought to be connected with pure CO or CO in a non- 
hydrogenated mixture containing species such as O2, N2 or 
CO2. The physical difference between the two components is 
that hydrogen-bearing species form pseudo-chemical bonds be- 
tween the hydrogen and the CO molecules, while the non- 
hydrogenated mixture only interacts via weak van der Waals 
forces. Traditionally, in the astronomical literature the two 
components are referred to as "polar" and "apolar". However, 
while it is true that the traditional mixtures used to repro- 
duce the red component are dominated by polar molecules, 
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Fig. 3. Sources in the Serpens cloud. Legend as in Fig.^ 



the chemical mechanism governing the profile of this com- 
ponent are bonds with the hydrogen atoms of the surround- 
ing H2O, CH3OH and NH3 molecules. Therefore, in the astro- 
nomically relevant case, the presence of hydrogenated species 
is what physically distinguishes the red component from the 



other two components and not their polarity. Incidentally, the 
presence of hydrogen bindings in the red component also deter- 
mines a number of known characteristics such as lower volatil- 
ity, ability to trap CO molecules in micropores etc. Henceforth 
in this article we will refer to the two components as hydrogen- 
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bonding and van der Waals-bonding components, which more 
accurately describes the chemical and physical structures in the 
ice mantles. 

A simple way of treating different components of the CO 
ice band, is to postulate that the entire band is fundamentally a 
superposition of contributions from CO stretching vibrations 
where the CO is situated on or in a discrete number of en- 
vironments. The term environment is defined in the broadest 
sense and refers to an average configuration of nearest molec- 
ular neighbours to a CO molecule for a given ice. The environ- 
ment depends on the abundance of other species as well as the 
structure of the ice, such as phase, density, porosity etc. Each 
environment should give rise to a certain spectral line with a 
center frequency and a width determined by the average inter- 
action between a CO molecule and its nearest neighbours. A 
detailed discussion and justification is given in Sec. [5] along 
with implications for the interpretation of the spectra presented 
in this paper. 

3.2. The analytical model 

To test the simple postulate on the large sample of YSO's ob- 
served here, every CO ice band is assumed to be decomposed 
into the three components mentioned in the literature. Since 
the presence of discrete components suggests that only a small 
number of environments are fundamental to ice on interstellar 
grains, it is also required that each of the three components has 
a constant central position and width such that only the rela- 
tive intensities are allowed to vary from source to source. In 
practice the position (but not the width) of the middle compo- 
nent is allowed to vary freely because it is the most sensitive to 



grain shape effects. The derived center position may then serve 
to test the degeneracy of the problem. It is essential to note, 
that if the fitted center positions of the middle component for 
the entire sample vary by more than is allowed by the statistical 
uncertainty of the data, then either the fit is degenerate or our 
assumptions are incorrect. 

Since there are sources in our sample where only one com- 
ponent in the CO ice dominates, it is straightforward to deter- 
mine the best centers and widths for each component. In this 
way TPSC 78, SVS 4-9, Cha INa 2 and RCrA IRS 7A are 
used to determine the position and width of the red compo- 
nent, while L1489, Reipurth 50, IRS 63 and Elias 32 fix the 
parameters of the blue component. 

The fundamental analytical shapes which can be chosen for 
each component are Gaussian and Lorentzian profiles. If three 
Lorentz curves are used in the fit, it quickly becomes clear that 
the blue side of the CO band fits very poorly in all cases. This 
is due to the steepness of the blue wing in the astronomical 
spectra. Since a Lorentz profile has a very broad wing, it is not 
appropriate for the middle and the blue component and conse- 
quently a Gaussian is adopted for these two components. At the 
same time, use of a Lorentz curve for the red component gives 
significantly better results than using a Gaussian for all three 
components. Consequently, two Gaussians and a Lorentzian 
were used for the blue, middle and red component, respectively. 
Physical reasons for this choice are given in detail in Sec. 14.21 

In addition to the three components used for the CO ice 
band, a Gaussian is placed at about 4.61 pm to account for pos- 
sible absorption by OCN . In this case our choice of a Gaussian 
is not as significant as for the blue and the middle component. 
A Gaussian gives in most cases a better fit than a Lorentzian, 
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Fig. 5. Sources in the Corona Australis cloud. Legend as in Fig.[0 



but due to the weakness of the feature other types of profile 
cannot be ruled out. This last component does not overlap with 
the CO band, and has rarely any influence on the three compo- 
nent fit and thus both center and width of this component are 
allowed to vary. 

Before attempting to fit the decomposition to the data a 
correction for the systemic velocity was applied to the model. 
The motion of the Earth in addition to the LSR velocity of the 
sources can shift the spectrum by up to +40 km s _1 correspond- 
ing to +0.29 cirT 1 . The fits improve significantly when a proper 
velocity correction is applied. Finally, the model spectra are 
convolved with a Gaussian to match the spectral resolution of 
the ISAAC spectra. Even at R = 5 000 (~ 0.4 cm" 1 ) this cor- 
rection must be done to obtain consistent results. 

The common presence of CO gas phase lines requires that 
care is taken when fitting a Gaussian to the 4.61 /mi band. 
Especially when broad emission lines are present, a careful se- 
lection of unaffected parts of the spectrum for the ice fit must 



be made. This is evident in e.g. IRS 43, Elias 32 and HH 100 
IR, and indeed in most of the sources to some extent. In gen- 
eral the CO gas-phase lines were identified and small areas be- 
tween the lines were picked by eye for the fitting. In Fig.[S]an 
example is shown of the line identification. This also prompts 
a concern in the case of spectra of low signal-to-noise, which 
may contain a significant filling-in of the 4.61 fim region from 
gas phase emission lines, which individually may remain un- 
detected. This may be the case in SVS 4-5, Cha IRN, Cha IRS 
6, RCrA IRS 7B and RCrA IRS 5B. 

Only the points on the spectrum which are not affected by 
CO gas-phase lines are used in the fit of the phenomenologi- 
cal ice band. In cases where broad emission lines with a com- 
plex structure are present, it is especially difficult to determine 
which points to use for the fit. This is in a few cases reflected 
in a reduced^ 2 value significantly greater than unity. However, 
since most of the broad gas lines are very weak compared to the 
main ice features, the final fitting parameters are not expected 
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Fig. 8. Example of the identification of gas phase CO emission 
lines. Shown is the spectrum of Elias 32. In this case all the 
lines are resolved. A weak and broad underlying 4.61 pm fea- 
ture is seen. 



to be severely affected. Note that the broad lines are fully re- 
solved and the selection of points will not improve with higher 
spectral resolution. 

3.3. The CO ice band structure 

The common parameters adopted for the decomposition are 
shown in Table|2] The best fit parameters are given in Tables[3] 
and @] and the decompositions are shown overlaid on the spec- 
tra in Figs. ^{7] The uncertainties on optical depths are given as 
statistical 3<x errors. Due to the systematic errors introduced by 
the continuum fit, the presence of gas-phase CO lines and the 
phenomenological decomposition we estimate that the statisti- 



cal uncertainties on the optical depths of the red, middle and 
blue components are underestimates of the total error when the 
features are strong. However, since the fitting parameters for 
the 4.61 pm feature are left free and the feature is not blended 
with other components, the statistical 1 <x error on r(2 175 cm -1 ) 
is probably a good estimate of the total error. For consistency 
reasons the errors on this quantity are also given as 3<x errors 
in Table E| 

The empirical fits are in general excellent, as is seen from 
the reduced x 1 values in Table and by visual inspection. 
The distribution of fitted positions of the middle component 
is shown in Fig. [H] It is found to vary between 2139.7 cm" 1 
and 2140.1 cm" 1 which corresponds to a standard deviation 
of 0.16 cm" 1 . Since the mean uncertainty of the positions is 
0.12 cm" 1 , the scatter in the fitted center positions is to a large 
extent caused by statistical uncertainties in the fitting and not 
necessarily by a real effect. An additional scatter may be in- 
troduced from uncertainties in the velocity correction of the 
model. The accuracy of the velocity correction is illustrated by 
the differences between the velocity shift of the CO gas-phase 
lines and the calculated velocity shift. If no velocity correc- 
tion is applied, the spread in center positions of the middle 
component increases significantly to a standard deviation of 
0.20 cm" 1 . The parameters for the blue and the red component 
are slightly more uncertain, since they are both shallower and 
the red one significantly broader than the middle component. 

In view of the large number of lines of sight probed, the 
large range of source parameters and the great variety in CO 
ice profiles, this result is highly surprising. In simple terms, it 
is seen that all CO bands towards both YSOs and background 
stars can be well fitted with an empirical model depending on 
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Fig. 7. Sources in Vela and the Taurus molecular clouds. Legend as in Fig. [2 



Table 2. Best fitting CO ice profile decomposition 



Component 


Center 




FWHM" 




[cm- 1 ] 




[cm- 1 ] 


red 


2136.5 


±0.3 


10.6 


middle 


2139.9 


+ 0.16 


3.5 


blue 


2143.7 


±0.3 


3.0 



' FWHM of Lorentzian profile for the red component and of 
Gaussian profile for the middle and blue components. 



only 3 linear parameters. Therefore our original postulate from 
Sec. B.ll that the CO stretching vibration mode can be decom- 
posed into a small number of unique components holds. This 
strongly suggests that there is more physics contained in the 
decomposition than previously assumed. The approach of fit- 
ting the same simple structure to every astronomical spectrum 
has thus provided three fundamental, identical line components 
that are common to all the lines-of-sight observed, and which 
must be taken into account in subsequent experimental and the- 
oretical modeling of interstellar ice mantles. 



OCN- ion JSchntte &Oreenber7lfl997tlPendleton et al]ll99'* 
iNovozamskv et alJbOOll) . The derived band parameters along 
with upper limits are given in Table |4] The widths and cen- 
ter positions of the band seem to vary significantly. In par- 
ticular, many of the bands are centered closer to 2175 cirr 1 
than to the ~ 2165 cirr 1 found for the bands discovered previ- 
ously, mostly towards high-mass sources. Also laboratory spec- 
tra give center positions of the OCN" band and for the CN- 
stretch, which in general are located red -ward of many of the 
obse rved center positions. For ex ample iHudson et alJ d200ll) 
and Ivan B roekhuiz en et alJ J2003h find that the centers of the 
CN-stretch in interstellar ice analogs vary between 2158 and 
2170 cm-'. Most of the fitted centers in our sample are placed 
between 2170 and 2180 cirr 1 . There is a clear tendency that the 
higher mass stars have redder 4.61 /im band centers than lower 
mass stars. The average position for the intermediate mass stars 
Reipurth 50, LLN 20, LLN 33 together with the low mass stars 
TPSC 1 and TPSC 78, which are located only 30" from the 
trapezium OB association in Orion, is 2167 + 2 cm-'. The rest 
of the detected bands have an average position of 2 1 73 +4 cirr 1 . 
For simplicity the band will hereafter be referred to as the 
2175 cm -1 band. 



3.4. The 4.61 /mi band 

A weak absorption band between 2200 cirr 1 (4.55 jjm) 
and 2150 cirr 1 (4.650 fim) is detected towards 30 of the 
39 sources. Absorption in this spectral region can be at- 
tributed to the stretching mode of CN bonds. The most likely 
species to contain the CN bond around young stars is the 



3.5. Correlations 

Correlation plots for the derived optical depths of the compo- 
nents of the decomposition are shown in Figs.l9UT3l Possible 
interpretations of the correlations, non-correlations and exclu- 
sion regions are given in Sec. 13.61 All correlations except Fig. 
^]only show points with firm detections of the relevant CO 
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Source 


Lada class 


Bolom. 


M-band 


Date of 


Standard star 


Int. time 


Simbad name 






lum. [L Q ] 


resolution 


observation 




[min] 




p Ophiuchus 


TO C /IT 

1Kb 42 


T TT 
1-11 


J.O 


c aaa 

3 uuu 


1 1 iq /oaa i 
13/0/2UU1 


rsSoUo4 


OA 

2U 


BKLl J loz /Z1-Z4414Z 


IRS 43 


1 


6.7 


10 000 


1 /C /OA AO 

1/5/2U02 


BS6175 


25 


ni/T rr T1£00"0£ 0/1/1AC1 

BKL1 J 162/26-244051 


IRS 44 


1 


8.7 a 


10 000 


1 /C /OA AO 

1/5/2UU2 


D C £ AO A 

BS6U84 


15 


1)1)1)/'/' A 1 ) 

BdKLu 49 


TD C A 

IKo 40 


1 


U.OZ 


1 n nnn 
1U UUU 


i /^ /o aao 
1/j/zUUz 


rSooUo4 


1 J 


T ) T ) I ) / ' / ' CO 


TD C A Q 

IKo 4o 




1 /111 

1 A 


1 n aaa 
1U UUU 


a /^ /onno 
4/j/zUUz 


rSooUo4 


OA 

2U 


r>KLl J loz lo /-z4jUjj 


tdc ci 
IKo 1 


1-11 


1 1 a 


^ nnn 

J UUU 


ly/o/zuui 


TDC^AQ/1 

rSoOUo4 


on 

ZU 


D17TT T 1 ^070A O A A O 1 f, 

DlvLl J 10Z / J7-Z44J 10 


IRS 54 




6.6 a 


10 000 


4/5/2002 


BS6175 


10 


BKLT J 16275 1-243 145 


IRS 63 


III 


4.2 b 


5 000 


19/8/2001 


BS6084 


20 


GWAYL 4 


WL6 




1.7" 


10 000 


4/5/2002 


BS6084 


30 


BBRCG 38 


WL 12 




2.6" 


10 000 


6/5/2002 


BS6084 


25 


BBRCG 4 


GSS30IRS1 




21 a 


5 000 


3/9/2001 


BS6084 


20 


BKLT J162621-242306 


CRBR 2422.8 


III 


0.36 a 


10 000 


5/5/2002 


BS6378 


35 


CRBR 2422.8-3423 


Elias 32 


II 


1.1* 


10 000 


1/5/2002 


BS6084 


25 


[B96] Oph B2 5 


VSSG 17 




3.7" 


10 000 


1/5/2002 


BS6084 


35 


[B96] Oph B2 6 


RNQ91 


II 


4.7 b 


5 000 


20/8/2001 


BS6147 


10 


HBC 650 


Serpens 


EC 90A 




35 d 


10 000 


1/5/2002 


BS6084 


16 


[B96] Serpens 7 


EC 90B 




20 d 


10 000 


1/5/2002 


BS6084 


16 


[B96] Serpens 7 


EC 82 


11 


9.7 d 


10 000 


1/5/2002 


BS6084 


25 


CK3 


SVS 4-9 




1.5 d 


10 000 


6/5/2002 


BS7121 


21 




SVS 4-5 




1.2 d 


10 000 


6/5/2002 


BS7121 


21 




CK2 


Background 


? 


10 000 


6/5/2002 


BS6378 


20 


CK2 


EC 123 


II 


l d 


10 000 


6/5/2002 


BS6378 


20 


[EC92] 123 


Chameleon 


ChaIRN 




15 e 


10 000 


7/5/2002 


BS3485 


30 


[AWW90] Ced 1 1 1 IRS 5 


Cha INa 2 




0.6 f 


10 000 


6/5/2002 


BS4844 


25 


[PMK99] IR Cha INa2 


Cha IRS 6A 




6 g 


10 000 


7/5/2002 


BS3468 


30 


[PCW91]Ced 110 IRS 6 A 


Corona Australis 


HH 100 IRS 




14 h 


10 000 


2/5/2002 


BS7121 


10 


V* V710CrA 


RCrA IRS5A 




0.6 e 


10 000 


2/5/2002 


BS7348 


20 


[B87] 7 


RCrA IRS5B 




0.4 e 


10 000 


2/5/2002 


BS7348 


20 


[B87] 7 


RCrA IRS7A 




3 1 


10 000 


2/4/2002 


BS6084 


30 




RCrA IRS7B 




3 1 


10 000 


2/4/2002 


BS6084 


30 




Orion 


Reipurth 50 




70) 


10 000 


12/11/2001 


BS3468 


20 


HBC 494 


TPSC 78 




~ 30 d 


10 000 


14/11/2001 


BS1790 


20 


TPSC 78 


TPSC 1 




~ 20 d 


10 000 


14/11/2001 


BS1790 


20 


TPSC 1 


Taurus 


LDN 1489 IRS 




3.7 k 


5 000 


19/9/2001 


BS1543 


20 


NAME LDN 1489 IR 


Vela 


LLN 17 




3100' 


800 


11/1/2001 


BS3185 


10 


IRAS 08448-4343 


LLN20 




344 1 


10 000 


14/11/2001 


BS3468 


8 


[LLN92] 20 


LLN 33 




91' 


10 000 


14/11/2001 


BS3468 


12 


[LLN92] 33 


LLN 39 




807' 


10 000 


14/11/2001 


BS1790 


20 


[LLN92] 39 


LLN 47 




21' 


10 000 


12/11/2001 


BS3468 


20 


[LLN92] 47 



Ref erence s teontemps eTail <200lh. b IChen et al ll995h . c IChen et alJ ll99%. d Rough estimate this wo rk.HChen et alj ll997l) . f IPersi et al 



1999h. 8 IPersi et alj <200ll) . " IWilking et alj 119*89^ ' IWilking et alJ <1986l) . > iDartois ;stsl\ j2003h . k iKenvon et aLr<1993h .' ILiseau et al 



1992) 



components. A decomposition of the CO ice band taken from 
ISO-SWS (A/AA ~ 2000) archival spectra of the high-mass 
stars W 33 A, GL 989, GL 2136, NGC7538 IRS 1 and RAFGL 
7009S have been added to the sample for comparison. The pa- 



rameters derived for these sources are included in Table fj and 
0] The error bars in the correlation plots are +l<r. 

No correlation is found between the red and the middle 
component (see Fig.|5}. Sources exist where virtually only the 
red component is visible, such as TPSC 78, TPSC 1 and Cha 
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Table 4. Best fits to the 2175 crrT 1 band 



Source 


center 




FWHM 


t(2175 cnr 1 )" 


p Ophiuchus 


IRS 42 


2175.7 


± 2.5 


23 ± 


6 


a a 1 o i a nnt 
0.012 ± O.UOo 


IRS 43 


2175.0 


± 1.2 


14 ± 


3 


A AT i A C\1 

0.07 ± 0.03 


IRS 44 


2174.8 


± 2.5 


11 ± 


6 


A A A i A AC 

0.04 ± 0.05 


IRS 46 


2165.0 


± 3.2 


24 ± 


8 


0.06 ± 0.03 


IRS 48 


2194.3 


± 8.6 


29 ± 


11 


AA11 i A A 1 C 

0.013 ± 0.015 


IRS 51 


2170.2 


± 1.2 


36 ± 


4 


A AOA i A f\f\C 

0.089 ± (J. 006 


TP Q *\A 










< U. 1 


IRS 63 


2176.2 


±0.8 


26 ± 


2 


0.073 ± 0.009 


WL6 










<0.06 


WL 12 


2175.2 


± 1.3 


11 ± 


3 


0.05 ± 0.04 


GSS30IRS1 










< 0.075 


CRBR 2422.8 


2175.2 


± 1.1 


9±2 


0.09 ± 0.05 


Elias 32 


2175.1 


±0.8 


29 ± 


2 


0.14 ±0.03 


VSSG 17 










< 0.2 


RN0 91 


2171.2 


±0.4 


28 ± 


1 


0.125 ±0.009 


Serpens 


EC 90A 


2170.5 


±0.7 


17 ±2 


0.073 ± 0.009 


EC90B 


2163.6 


±0.4 


10 ± 1 


0.10±0.02 


EC 82 










<0.06 


SVS 4-9 


2171.3 ± 1.9 


14 ± 


4 


0.15 ±0.10 


SVS 4-5 


2166.6 ±3.4 


19 ± 


8 


0.10±0.10 


CK2 










<0.15 


EC 123 










<0.3 


Chameleon 


ChaIRN 


2168.6 


± 3.0 


22 ± 


6 


0.13 ±0.07 


Cha INa 2 










<0.09 


Cha IRS 6A 


2169.1 


±2.1 


14 ± 


4 


0.15 ±0.14 


Corona Australis 


HH 100 IRS 


2173.2 


±0.2 


19 ± 


1 


0.072 ± 0.006 


RCrA IRS5A 


2172.7 


±0.9 


16 ± 


2 


0.08 ± 0.02 


RCrA IRS5B 










<0.09 


RCrA IRS7A 


2171.5 


±2.2 


23 ± 


5 


0.10 ±0.06 


RCrA IRS7B 


2167.0 


± 1.2 


12 ± 


2 


0.15 ±0.08 


Orion 


Reipurth 50 


2170.2 


±0.1 


27 ± 


1 


0.129 ±0.003 


TPSC 78 


2165.3 ±0.2 


20 ± 


1 


0.111 ±0.006 


TPSC 1 


2167.9 ±0.7 


26 ± 


2 


0.18 ±0.03 


Taurus 


LDN 1489 IRS 


2176.9 


± 1.2 


29 ± 


3 


0.09 ± 0.02 


Vela 


LLN 17 


2169.2 


±0.3 


27 ± 


1 


0.13 ±0.01 


LLN20 


2166.2 


±2.6 


31 ± 


9 


0.09 ± 0.04 


LLN 33 


2167.4 


±0.7 


17 ± 


2 


0.19 ±0.04 


LLN 39 


2173.4 


± 1.3 


24 ± 


3 


0.010 ±0.003 


LLN 47 










<0.03 


Additional sources 


GL 2136" 


2164.4 


±0.2 


24 ± 


2 


0.13 ±0.01 


NGC 7538 IRS 1* 










<0.03 


RAFGL 7009S'' 


2168.2 


±0.2 


25 ± 


1 


0.47 ± 0.02 


W33A* 


2165.7 


±0.1 


26 ± 


1 


1.19 ±0.03 


RAFGL 989* 


2169.7 


± 1.2 


23 ± 


3 


0.04 ± 0.01 



" 3tr errors are given. 
b ISO-SWS archive data. 




CK 2 



SVS 4-9 

I ^ I 

RAFGL7009S. 



i.O 0.2 0.4 0.6 0.8 1.0 1.2 

CO 2136.5 cm -1 (re) optical depth 

Fig. 9. Red component optical depth against middle component 
optical depth. The solid line corresponds to 16% of the total 
amount of CO present in a hydrogen-bonding environment un- 
der the assumption that the red component is carried by a water- 
rich mixture. The percentage can be calculated using Eqs.[3]and 

El 



INa 2, while other sources are dominated by the middle com- 
ponent. However, no source with a deep middle component has 
been found where the red component is altogether absent. This 
effect is seen in the figure as an exclusion region in the upper 
left corner where no points are located. The region is indicated 
by a solid line. 

In Fig.llOlthe red component optical depth is plotted against 
the ratio of the middle and red component optical depths. 
Evidently no direct correlation is seen. However, the region in 
the upper right corner of the figure is excluded, which means 
that a deep red component results in a low strength of the mid- 
dle component relative to the red component. This plot has 
been included to investigate if the ratio between the two com- 
ponents can be used as a processing indicator. 

The plot of the red and the blue component shown in Fig. 
ITT1 also shows that certain regions in parameter space are ex- 
cluded as indicated by solid lines. However, it should be noted 
that the blue component is the most susceptible to systematic 
effects created by the artificial decomposition due to the domi- 
nating presence of the nearby middle component. Thus caution 
should be exercised when interpreting this correlation plot in 
areas where the blue component is weak. In essence, a strong 
blue component is excluded if the red component is weak as 
seen by the absence of sources in the upper left corner of the 
figure. The absence of sources in the lower right corner indi- 
cates that a deep red component excludes a weak blue compo- 
nent 

Fig. 1121 shows the depth of the middle component plotted 
against the blue component. The curves indicate constant po- 
larisation fractions of the background source as required if the 
blue component is assumed to be due to the LO-TO splitting of 
crystalline CO as discussed in Sec. 14.71 

Finally, a trend between the red component and the 
2175 cm" 1 component is evident in Fig.^] In this correlation 
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Fig. 10. Red component optical depth against the ratio of the 
middle component optical depth and the red component optical 
depth. The solid arrow indicates the possible evolutionary track 
of a source if the pure CO does not migrate into a porous water 
ice during warm-up. The solid line indicates the excluded re- 
gion. The dashed curves indicate evolutionary tracks for SVS 
4-9 assuming that migration does occur with efficiencies of 0. 1 
and 0.5. The dot-dashed curves are evolutionary tracks for WL 
6 with migration (see Sec. 13. 61 ). 




CO 2136.5 cm 1 (rc) optical depth 

Fig. 11. Red component optical depth against blue component 
optical depth. The solid lines indicate excluded regions. 

plot the sources have been divided into low-mass and high- 
mass stars to search for evidence for enhanced OCN abun- 
dances around high-mass stars. High-mass stars have been de- 
fined as stars more luminous than 50 L (bolometric). This def- 
inition is flawed by the fact that the luminosity of a young star 
is dependent on age as well as inclination, but due to the in- 
herent difficulties in determining spectral classes of embedded 
stars no better definition was found. Roughly, the definition will 
divide bona fide embedded sources into young stars lighter and 
heavier than ~ 2 M , respectively. The low-mass stars exhibit a 
nearly linear relation between the red component optical depth 
and the 2175 crrr 1 feature for r(rc) up to 0.6 after which the 
relation apparently saturates and flattens out. However, the flat- 
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Fig. 12. Middle component optical depth against blue compo- 
nent optical depth. The solid lines are curves of constant po- 
larisation fraction calculated using Eq. [5] The blue compo- 
nent optical depths have been slightly corrected by using a 
Lorentz oscillator model for the middle component rather than 
a Gaussian. See Secs. l4.2l and l4.7l 
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CO 2136.5 cm 1 optical depth 

Fig. 13. Optical depth at 2175 cm -1 against optical depth of the 
red component. Stars indicate low-mass young stars defined by 
having Lboi < 50 L . Diamonds indicate intermediate to high- 
mass stars with Lboi > 50 L . Small arrows indicate upper lim- 
its on the 2175 cirT 1 optical depth. The points of the high-mass 
stars W 33A and RAFGL 7009 S lie outside the plotting range 
and are indicated by large arrows. Although TPSC 1 and TPSC 
78 are low-mass stars, they are marked as high-mass stars due 
to their location in the high-mass star forming Orion core. The 
marked region indicates the possible range of the correlation 
for low mass stars. 

tening in the relation only depends on a few points and should 
thus be interpreted with caution. Some of the high-mass stars 
follow the same relation, but many show a dramatic excess ab- 
sorption around 2 1 75 cirT 1 . When including the massive young 
stars W 33A and RAFGL 7009S a rough correlation can also 
be seen for the high-mass stars, although much steeper than 
for the low-mass stars. No sources are found significantly un- 
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der the best-fitting line and no 2175 cirT 1 absorption bands are 
detected for sources without CO ice bands. Within our detec- 
tion limits, there seems to be an almost one-to-one relationship 
between the presence of a 2175 cirr 1 feature and the red com- 
ponent of the CO ice. 



3.6. Possible interpretation of correlations 

In Sec. 13.51 the optical depths of the four separate solid state 
components identified in the spectra have been plotted against 
each other. The regions of the parameter space covered by the 
observed sources and especially the regions where no sources 
fall may provide clues to the nature and evolution of the com- 
ponents. 

If the middle component is ascribed to pure CO and the 
red component to CO in a water-rich environment, one rele- 
vant question is how the volatile, pure CO interacts with the 
water-rich CO when the ice is thermally processed. If pro- 
cessing causes the evaporation of some of the carrier of the 
middle component while the rest migrates into a porous water 
ice belie ved to be located below, such as has been suggested 
bv lCollings et alJ J2003). an increase in red component optical 
depth would be expected when the ratio of middle to red com- 
ponent decreases. Due to different absolute column densities 
along different lines of sight, this effect may be hard to detect 
without comparing to an independent ice temperature indicator, 
such as the 3.08 fim water band. However, the exclusion regions 
seen in Figs.[9]and^)|still provide useful constraints. The fact 
that the sources with the deepest red components seem to have 
the smallest middle components compared to the red as shown 
with a solid line in Fig.^|indicates that this redistribution of 
CO may take place. If this holds then Fig.|9]would indicate that 
some redistribution of CO has taken place for all sources, since 
the red component is always present with a minimum depth 
relative to the middle component. This is indicated in Fig. [9] 
with a solid line, corresponding to 16% of the total amount of 
CO embedded in water. The interpretation is then that the least 
processed CO profiles are those closest to the solid line, and for 
some reason the minimum fraction of CO in water along any 
line of sight is 16%. 

Sketches of possible evolutionary tracks with or without 
migration of the CO ice into the water ice during warm-up 
have been drawn in Fig. ^3 The solid track shows the mid- 
dle component evaporating independent of the red component 
at temperatures below 90 K for SVS 4-9. This scenario of no 
CO migration raises the question of why no "progenitors" have 
been found of the sources in the lower right corner of Fig [K)] 
such as SVS 4-9, RCrA IRS 7B and the massive YSOs W 33A 
and RAFGL 7009S. According to the simple evolutionary track 
without migration, these progenitors should be found in the up- 
per right corner of the plot. If migration does occur, i.e. if the 
CO of the middle component is deposited on a layer of porous 
water ice, the first part of the evolutionary track changes due to 
a growth of the red component simultaneously with the evap- 
oration of the middle component. For a constant efficiency, E, 



of migration, such that dN rc = —EdN mc , the evolutionary track 
is: 

Nmc_ _ 1 Wc,0 - Nrc) + N mc fl 

N rc ~ E N rc 

where N mc and N IC are the column densities of the middle 
and red components, respectively, while N mc fl and N IC ,o are the 
observed column densities. In Fig. EH evolutionary tracks for 
migration efficiencies of 0.1 and 0.5 have been drawn through 
the extreme sources WL 6 and SVS 4-9. It is seen that the CO 
profile of SVS 4-9 can be explained by evolution from a quies- 
cent line of sight, such as that probed by the background source 
CK 2, if the migration efficiency is large. This is consistent 
with the lack of sources above the solid line in the diagram. 
The tracks through WL 6 show how this source can evolve into 
sources in the lower left part of the diagram. The distribution 
of sources in Fig.[TOlcan thus be taken as evidence that migra- 
tion of CO molecules into a porous water ice does occur upon 
warm-up of the circumstellar grains (but see Sec. l4.5> . 

Amorphous water ice only becomes porous when it is de- 
posited at low temperatures (10-20 K). The migration effi- 
ciency drops with increasing deposition temperature and if the 
ice is deposited at temperatures higher than 70 K, the pores 
do not form at all JCollings et al.lEo03l) . If the distribution of 
sources in Fig.[H)]is interpreted as a result of migration, an ice 
structure created at a low deposition temperature is favoured. 
Since water is thought to be formed on interstellar grain man- 
tles rather than being deposited, this may put constraints on the 
formation of H2O, since a similar porous structure must then 
be formed as the water is chemically assembled on the grain 
surfaces. 

The trend between the red component and the 2175 cm 
band has not previously been observed, probably because ear- 
lier samples contained a large fraction of intermediate and 
high mass stars. The observation that high-mass stars exhibit 
a radically different correlation suggests that two different 
components contribute to the absorption between 2155 and 
2185 cur 1 . Since the 2175 cirr 1 profiles towards the low mass 
stars also tend to be centered redder than laboratory experi- 
ments allow for the CN-stretch, we propose that a second weak 
absorber is present with a center at 2170 - 2180 cirr 1 and a 
FWHM of 15 - 25 cirr 1 and that, based on the correlation with 
the red component, a likely carrier is CO in a presently uniden- 
tified binding site. One significant implication of the proposed 
correlation is that no absorption at 2136.5 cirr 1 is expected 
for any source if the 2175 cirr 1 band is not present. This is 
supported by the sensitive upper limits given in Fig. [O]for a 
number of sources which contain abundant water ice. 

The excess seen in some high mass stars is then an indica- 
tion of the formation of OCN by thermal, UV or other forms 
of processing not present in low mass stars. The OCN band in 
these sources floods the weaker underlying band and the corre- 
lation disappears. 

A number of alternative explanations may be explored. For 
example, it is known that a libration band of water ice creates a 
broad band centered on 2200 cm . This libration mode could 
create an excess absorption, which when requiring the contin- 
uum to fit the blue edge of the spectra, may mimic the observed 
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Fig. 14. Distribution of fitted middle component positions. The 
thick, dashed curve shows the distribution derived from the 
astronomical spectra and corrected for the systemic veloci- 
ties of the sources. Only middle component positions with lcr 
errors of less than 0.15 cm 1 have been included. The thin 
curve shows the distribution of the center positions of the non- 
hydrogenated ices from the ice analog database by E97. Only 
ice analogs with FWHM between 3 and 4 crrT 1 have been in- 
cluded in the distribution. All the laboratory ices are corrected 
for grain shapes as indicated in the figure. 

2175 cirT 1 band. This is a valid objection, since the depth of 
the feature depends on the assumption that the blue edge of the 
spectrum can be used as continuum. However, the water libra- 
tion band has a FWHM of about 300 crrT 1 and the total depth is 
only ~ 2.5% of the 3.08 /mi main water band. This means that 
the expected excess optical depth from the water combination 
mode when using a continuum fixed at 2100 and 2200 cirT 1 is 
at most t ~ 0.01 for a r ~ 5 main water band. 

4. Physical modeling of the solid CO band 

4.1. The middle component 

Having identified a simple and well-defined decomposition of 
the astronomical CO ice band, a unique basis is provided, 
which can be tested against models of interstellar CO ice. Since 
the shape of the middle component is the best constrained, 
it is instructive to compare it to laboratory spectra of van- 
der-Waals interacting ice mixtures corrected for different grain 
shapes. Such a comparison is not intended to provide unique 
constraints on the composition of the CO ice responsible for 
the middle component, but illustrates well how narrow the ob- 
served range of profiles is, compared to the possible range of 
shapes depending on ice composition and grain shape correc- 
tions. For this purpose the optical constants from the laboratory 
databa se of non-hydrogenated mixtures by lEhrenfreund et alJ 
il 19971) (hereafter E97) are used. The laboratory spectra are 
corrected for grain shape effects using the standard four par- 
ticle shapes almost exclusively used in the literature, namely 
spheres, coated spheres, a continuous distribution of ellip- 
soids (CDE) and an MRN jMathis etal.lll977h size distribu- 
tion of coated spheres. The parameters for the grain shape cor- 



rection from E97 have been adopted. For details the reade r 
is referred to iTielens et alJ Jl99ll) : lEhrenfreund et alJ dl997l) : 
iBohren & Huffmanld 19831 In general, the grain shape correc- 
tions are independent of grain size as long as the grains are 
small, i.e. 2na/A <K 1, and the Rayleigh approximation is 
valid. Clearly, the requirement is only marginally satisfied at 
5 fim and breaks down for grains larger than 0. 1 fim. This may 
be a serious concern in the vicinity of protostars, where grain 
growth is expected to occur. Strictly, grain shape corrections 
for wavelengths shortwards of 10 fim should be applied using 
a full numerical solution of the field interactions with the parti- 
cle. Nonetheless, we assume in the following that the Rayleigh 
condition holds given the limited constraints on the detailed 
grain size distributions in the observed lines of sight. In essence 
it is assumed that the extinction is dominated by small grains. 
This also implies that the effects of scattering out of the line of 
sight are ignored. For small particles, the ratio of average ab- 
sorption to scattering cross sections can be approximated with 
~ (18^/3) x v~ 3 V~', where v is the wave number and V is 
the particle volume, which only becomes less than unity for 
a > 10 fim at a frequency of 2140 cm -1 . For the CO band 
it is a good approximation to ignore scattering effects, while 
they become important for bands at higher frequencies, such as 
for the water band at 3 fim. Scattering effects on the CO band 
may become observable in regions with grain growth, such as 
in circumstellar disks. 

The coated sphere model was calculated using identical 
volumes of core and mantle, while the MRN model was cal- 
culated using a constant mantle thickness of 0.01 fim; mod- 
els with thicker mantles approximate the results for solid ice 
spheres. For these models, a constant baseline was subtracted 
from the absorption cross sections to remove the absorption 
from the silicate core. 

A Gaussian was fitted to the CO band to determine the peak 
position and FWHM of the laboratory ice absorption. Shape 
corrected laboratory ices with FWHM larger than 4 crrT 1 and 
smaller than 3 cirT 1 were discarded since they are inconsis- 
tent with the astronomical spectra. The distribution of the fitted 
center positions is compared to the distribution of center po- 
sitions of the laboratory spectra in Fig.^J Only five different 
laboratory spectra of non-hydrogenated ice mixtures are in the 
narrow parameter range defined by the observed middle com- 
ponent using the standard grain shapes out of a possible 228. 

The variety of laboratory mixtures that are consistent with 
the observational peaks range from pure CO at 10 K with a 
CDE grain shape, a H 2 0:CO:0 2 = 1:20:60 at 30 K for coated 
spheres, to H 2 O:CO:O2:N 2 :CO2=l:50:35:15:3 at 10 K for all 
the grain shapes. The latter mixture gives a reasonable center 
and width for all grain shapes since the dielectric function is 
weak. The shift and the broadening are therefore caused by 
the presence in the mixture of C0 2 and 2 , respectively, rather 
than by grain shape effects. It will be shown in Sec. l4.4l that the 
shape of the 13 CO band rules out the multi-component mix- 
tures. This range in possible mixtures illustrates the degener- 
acy inherent in fits to the middle component. However, it also 
shows that the range of ice mixtures in which the carrier of the 
middle component is present must be very limited since most 
of the non-hydrogenated ice analogs are excluded, i.e. the lab- 
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oratory spectra coupled with commonly used grain shape cor- 
rections show a much larger range in centers and band widths 
than those commensurate with the observations. This is encour- 
aging, since it points towards a simple solution to the questions 
regarding the composition and structure of CO-rich ice on in- 
terstellar grain mantles. 

We will proceed to show that the simplest physical model of 
pure CO ice is sufficient to explain all the middle components 
observed. Furthermore, it will be shown that in nearly all cases 
perfect fits are achieved using the CDE grain shape. 

4.2. Lorentz oscillators as models for binding sites 

The simplest possible physical model of the absorption of a 
molecule in the solid phase is that of a Lorentz oscillator. This 
is fundamentally a classical model but yields expressions sim- 
ilar to those of simple qua ntum mechanical models (see e.g. 
lGadzukll987tlZimarJ 19791) . albeit with significant conceptual 
differences. Interactions with the surrounding medium may of 
course alter the potential energy surface of the oscillator, af- 
fecting the line shape. 

In the classical picture the CO molecules are seen as a 
set of identical springs with mass m, corresponding to the re- 
duced mass of the CO molecule, charge e, spring constant K 
and damping constant b. The last two parameters are affected 
by the bond of a CO molecule to the surrounding molecules, 
and for each configuration of a CO molecule and its nearest 
neighbours a unique set of K and b exists. In an amorphous 
ice there is a continuum of different configurations for a clas- 
sical CO molecule, which when averaged over a large num- 
ber of CO molecules results in a characteristic set of K and b. 
In the simplest model the configurations of the CO molecules 
are random, resulting in a "pure de-phasing" broadening. See 
Appendix|X]for the technical argument for modeling solid state 
dielectric functions with Lorentz Oscillators. 

The complex dielectric function for a Lorentz oscillator is: 

e = eo + — £— , (2) 

0)^ - u) z — iyu 

where ix? p - e 2 N/mso is the plasma frequency, ot^ = K/m 
and y = b/m. N is the number density of oscillators. Finally 
eo determines the dielectric function at frequencies which are 
low compared to the electronic excitation frequencies. The pa- 
rameter is basically the low-frequency wing of the combined 
dielectric functions of all the electronic transitions. 

A thick slab, i.e. much thicker than a monolayer, of pure 
amorphous CO ice is expected to show just a single dominant 
environment, namely CO on CO, and the simple model should 
fit well to the measured optical constants. In the case of CO 
this is well known: the vibrational spectrum of adsorbed CO 
is commonly us ed as a probe o f surfaces in the chemical lit- 
erature (see e.g. Somorjai 1994, and references therein). Sub- 
monolayer coverages of CO can be used to probe the underly- 
ing surface structure, physical and chemical behaviour of CO 
on the surface, interactions between CO and other adsorbates at 
the surface, and to identify the range of binding sites favorable 
to CO on any particular surface. 



Note that dielectric functions superpose such that addi- 
tional binding sites and environments can be included in the 
model simply by adding single Lorentz oscillators in e-space. 
In this picture, the optical constants of more complicated ice 
mixtures can also be reproduced by fitting a sufficient number 
of Lorentz oscillators accounting for specific environments in- 
cluding other ice species, although the problem may quickly 
become too degenerate to yield useful physical information. 

Naturally, this is a very simplified picture and many other 
effects may play a role. For instance, some hydrogenated 
molecules, such as H2O, complicate the picture significantly by 
forming hydrogen bonds with the CO molecules, thus breaking 
the assumption of pure de-phasing. However, for weakly inter- 
acting environments such as that of CO interacting with itself 
or with N2 or O2, such a model may yield useful information. 
This is a break with the traditional view that solid state features 
from amorphous ices have a continuous range of shapes de- 
pending on the abundance of other species, but should be easy 
to verify experimentally using high resolution laboratory spec- 
troscopy in the limit of a low concentration of contaminating 
molecules, which interact with CO molecules only through van 
der Waals bonds, such as O2, N2 and CO2. In the following 
the term "environment" means a sum of unique sets of nearest 
neighbours giving rise to a sum of Lorentz oscillators. We will 
show that in the case of pure CO a simplified physical model 
can be used to fit the data. 

4.3. Oscillator density as fitting parameter 

There is a standing controversy in the literature regarding the 
accuracy of the determination of optical constants from labo- 
ratory measurements. Since the grain shape correction is sensi- 
tive to small differences in optical constants, the derived band 
profiles of mixtures with a high concentration of CO can vary 
significantly when using differe nt optical co nstants from the 
literature (e.g. E97 and Boog ert et al]l2002al) . This is partic- 
ularly true for pure CO. The resulting differences in the di- 
electric functions are illustr ated in Fig . H31 whe re the optical 
constants for pure CO f rom Baratta & Palumbo (1998) (here- 
after BP) and lElsila et all dl997h Thereafter EAS1 have been 
used to calculate the dielectric functions in panels (b) and (c), 
respectively. The main difference between the functions is in 
the plasma frequency, resulting in strengths of e varying by 
more than 20%. Also, the fits to a single Lorentz oscillator are 
slightly worse for the dielectric functions from BP and EAS. 
There is a second component clearly present in the spectrum 
from EAS with coq = 2142.3 cm" , which has been included 
in the fit. Thi s may be a con taminant. Also, recent results and 
work by ICollines et al.l d2003l) clearly shows the presence of 
this second peak, attributed to LO-TO splitting. Further dis- 
cussion of this effect is found in Sec. 14.71 of this article. The 
width (y) and center positions {ojq) of the dielectric functions 
are very similar in all the experiments. Within the model, the 
only parameter which can change the plasma frequency, and 
thus the optical constants, is the number of oscillators per unit 
volume, N. In other words, the difference between the labora- 
tory spectra of identical mixtures can to a large extent be ex- 
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Fig. 15 . Fit of a single Lorentz os cillator to the dielectric function of a pure CO ice at 10 K from E97 (al. iBaratta & Palumbo 
( 1998) (b) and lElsila et all i 19971) (c). The best fitting parameters are: (a): y — 1.5 cm , a> p — 195 cm 1 and eo = 1.67. (b): 
y = 1.5 crrT 1 , a> p = 160 cm -1 and eo = 1-57. (c): y = 1.75 cirT 1 , u> p = 175 cm" 1 and eo = 1-67. All have a>o = 2138.5 cm" 1 
Additionally a second component has been fitted in panel (c) with parameters: y = 1.2 cm" 1 , a> p = 30 cm" 1 and ojq = 2142.3. 



plained by differing ice density, differing ice porosity, or, less 
likely, by dilution by a relatively inactive molecule such as N2. 
Uncertainties in the determination of ice thickness in lab ex- 
periments may also play a role in the differences of the optical 
constants obtained by different groups. 

Fits using a single Lorentz oscillator as a model of the mid- 
dle component to the interstellar spectra with the highest S/N 
ratios are shown in Figs. 1161 and!17l The model uses the same 
phenomenological components for the red and blue compo- 
nents, but now adopts the Lorentz oscillator fit to the pure CO 
spectrum by E97. The middle component is grain-shape cor- 
rected using CDE particles. The CDE particles are used sim- 
ply because they fit the data extremely well and because this 
model is believed to simulate irregular particles well by in- 
cluding a maximal range of possible shapes. The CDE parti- 
cles have a specific distribution and have thus no free parame- 
ters. The actual shapes of interstellar grains are expected to be 
entirely different. The plasma frequency is allowed to vary to 
account for different ice densities or for small differences in di- 
lution. u) p was found to vary slightly between values of 170 to 
180 cm" 1 corresponding to densities 8 to 13% smaller than the 
E97 laboratory ice but very similar to the densities of the BP 
and EAS ices. Overall the fits are excellent and often clearly an 
improvement over the Gaussian model. Note that the fact that 
the grain-shape corrected Lorentz oscillators so closely emu- 
lates Gaussian profiles indicates that the blue and middle com- 
ponents are ices with a high concentration of CO, thus making 



grain shape effects important, while the red component is likely 
much more dilute, since it is close to the Lorentzian profile. 

It is instructive to repeat the exercise for the other three 
grain shapes used in the literature. Fig. ^] shows the CO pro- 
file for IRS 43 compared with the pure CO profile corrected 
for spherical grains, coated spheres and coated spheres with an 
MRN size distribution. Clearly none of the other Rayleigh limit 
grain shapes fit the observed profile. If any of these grain shapes 
is to be made to fit, some ice mixture is needed to broaden the 
profile. Alternatively, the MRN size distribution can be modi- 
fied. Since the MRN distribution is a power law with exponent 
-3.5, it is dominated by small grains. This is reflected in the 
CO profile by the fact that it peaks at 2140.2 cm" 1 ; i.e. a large 
mantle to core volume ratio shifts the profile to the blue. If a 
power law grain size distribution is to fit the red wing of the ob- 
served profile, the distribution must be much shallower in order 
to include a larger fraction of large grains. A broader shape- 
corrected profile can be obtained if the exponent is changed to 
a value closer to -1. This may be interpreted as evidence for 
grain growth near the young stars. However, this scenario re- 
quires that a change in profile is seen in YSO's as compared to 
quiescent dense clouds. Fig. ^] shows the observed CO pro- 
file of the well-studied background star CK 2. This is most 
likely a K0/G8 supergiant seen thro ugh a dense part of the 
Serpens molecular cloud core JCasali & Eiroall9 96). While the 
CO profile towards CK2 is severely affected by absorption lines 
intrinsic to the star as well as being saturated, it is clear that the 
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Fig. 16. Fit of a pure CO Lorentz oscillator model corrected for CDE grain-shape effects to the best observed spectra which are 
dominated by the middle component. The red and blue components have been added to the shape-corrected Lorentz oscillator 
as in the purely phenomenological model (cf. Table [3}. Often, the depth of the blue component had to be lowered by up to 
30% for sources with a deep middle component when using a Lorentz oscillator due to the presence of a blue wing in the CDE 
corrected profile. For the parameters of the adopted Lorentz oscillator the fit to the optical constants from E97 are used, except 
for the plasma frequency which is allowed to vary from 170 - 180 cm -1 . The dashed curve shows the sum of the red and blue 
components. A dashed straight line indicates the continuum level at zero optical depth. 
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Fig. 17. Same as Fig.^]for sources dominated by the red and blue components. For these sources a clear deviation from the pure 
CO model is seen between 2140 cirT 1 and 2143 cm -1 . 




Fig. 18. Comparison of different grain shape corrections for 
a Lorentz Oscillator pure CO. Solid line: CDE; Dotted line: 
Identical spheres; Dash-dotted line: identical coated spheres 
with equal volume core and mantle; Dashed line: MRN dis- 
tribution of coated spheres with a 0.01 ^m mantle thickness. 



wings of the CO profile are consistent with the shape given by 
CDE grains and inconsistent with MRN grains. In conclusion, 
the data presented here do not provide evidence for a strong 
difference in grain shapes in the quiescent medium compared 
to lines of sight towards YSOs. The similarity between the CO 
profile towards CK 2 and the other observed profiles also show 
that the separation of a possible contribution from solid CO in 
quiescent foreground clouds may be very difficult using only 
the solid state band. Gas-phase observations of rotational lines 
are thus needed to estimate the foreground contribution to a 
CO ba nd towards a YSO in a rigorous manner jBoogert et alJ 
l2000bl) . 

The excellent fit to the middle component as a single CO 
environment is strong evidence that the non-hydrogenated or 
van-der-Waals interacting component of interstellar ices is pure 
CO with at most a 5-10% contaminants. This percentage cor- 
responds to the range of oscillator strengths of the pure CO 
dielectric function (from 170 to 180 cirT 1 ) used to fit the data. 
With the current lack of knowledge about the true grain shapes 
in the interstellar medium, it can be concluded that within the 
quality of the presented data, the middle component is indistin- 
guishable from the profile of pure CO along all observed lines 



Fig. 19. Model profiles using a Lorentz oscillator for pure CO 
at 10 K corrected for CDE grains (solid curve) and MRN- 
distributed coated spheres (dot-dashed curve) compared to the 
spectrum of CK 2. 

of sight where the middle component dominates the CO profile. 
Deviations are seen in sources dominated by the red compo- 
nent. This is most clearly evident in Reipurth 50 and L 1489 as 
shown in Fig.[n] While the pure CO component is still present, 
one or perhaps two new features appear around 2141 crrT 1 and 
2142 cirT 1 . The features may also be present in IRS 63. These 
features have optical depths of only 0.1-0.2 and may therefore 
be present in all spectra, but are simply swamped by the deep 
middle component when this dominates. The centers of these 
new components coincide well with the bulk of the laboratory 
spectra in Fig.^] and may be evidence for a small contribution 
from mixed CO ices. 

For the pure CO combined with CDE grains the column 
densities can be found from: 

Npure.CDE = 6.03 CWT 1 X T max X A^, (3) 

where Ab u ik is the band strength of the bulk material. The nu- 
merical factor takes into account the r = 1 band equivalent 
width and the fact that the grain shape effect changes the effec- 
tive band strength with respect to the band strength measured 
in the laboratory with a factor of 0.71 = J Ccde^v/ J C^dv 
in this specific case. Ccde and C s i a b are the absorption coeffi- 
cients for pure CO CDE grains and for a pure CO slab, respec- 
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Fig. 20. The solid 13 CO band of IRS 51. The upper solid line 
is a Lorentz oscillator fit for a pure CO ice. The blue wing 
is somewhat affected by a broad CO gas phase emission line. 
The middle and lower curves are laboratory spectra of binary 
mixtures at 10 K rich (1:1) in O2 and N2, respectively. The 
vertical line indicates the center of the pure CO profile. 

tively. In the following we use a labo ratory band strength of 
1.1 x 10" 17 cm molec" 1 JGerakines et alll995l) . 

4.4. 13 CO in IRS 51 

Solid I3 CO is detected towards one of the sources, namely IRS 
51 (Fig. EH). Since the isotopic ratio of 12 CO to 13 CO is ex- 
pected to be between 65 an d 75 in the gas phase fo r typical 
molecular cloud conditions tLanger ^Pgnzia jll993h a nd be- 
tween 55 and 85 in the solid phase JBoogert et al.l2000ah . 13 CO 
is highly dilute in the ice. An important consequence is that 
grain shape effects are unimportant and the solid I3 CO band 
therefore offers the chance to disentangle, with respect to the 
band p rofile, the effects of th e ice matrix from those of the grain 
shape JBoogert et al.l2002al) . 

For the 13 CO band towards IRS 51 the same procedure 
as for the middle component of the I2 CO band is repeated. 
A Lorentz oscillator was fitted to the dielectric function by 
E97. The best fitting parameters are u> p = 22 crrT 1 , wo = 
2092.3 crrT 1 and y — 1.5 cm -1 . Again the plasma frequency 
is adjusted to fit the observed band. A good fit to the inter- 
stellar profile is found with lo p — 21 cm -1 , although the blue 
wing contains flux from the v = 2 — 1 P(6) (2193.40 cm' 1 ) 
and the v = 1 — P(12) (2194.86 crrT 1 ) emission lines from 
gas phase I2 CO. The red wing of the observed band is fully 
consistent with a pure CO ice. The ratio of the integrated 
imaginary part of the Lorentz oscillator dielectric functions 
gives the isotopic ratio of solid CO along the line of sight, 
AT( 12 CO)/jV( 13 CO) = 68 + 10 for IRS 51. The uncertainty re- 



flects the fitting uncertainty in the saturated I2 CO band assum- 
ing that the shape of the band is identical to the other sources. 
Comparisons with laboratory spectra of binary mixtures with 
O2 and N2 at 10 K from E97 are also shown in Fig.|20] Both 
laboratory mixtures give significantly worse fits to the IRS 5 1 
spectrum, the 02-rich mixture being too broad, and the N2-rich 
mixture being shifted to the blue. Heating the 02-rich mixture 
to 30 K narrows the profile, but not enough to equal the quality 
of the fit of that of pure CO. It may also be considered unlikely 
to find most of the volatile solid CO component in dense clouds 
at a temperatures higher than 20 K. 

The results for th e solid 13 CO band agre e well with the con- 
clusions reached by Boogert et al. (2002a) for the I3 CO band 
towards the high mass source NGC 7538 IRS9, who also found 
both isotopic bands of CO to be fitted well with pure CO. 
Furthermore, the isotopic ratio for NGC 7538 IRS9 of 7 1 + 1 5 is 
remarkably similar to that found for IRS 5 1 . Indeed, the shapes 
of the two 13 CO bands match closely, in agreement with the 
similarities observed among the I2 CO bands. 

4.5. The red component 

The red compone nt in the i nterstellar CO ice profile was fi rst 
assigned by Sandford et al. ( 1988) and Tielens et al. ( 199 JJ) to 
CO mixed in a hydrogen-bonding "polar" mixture containing 
species such as H2O or CH3OH. The previous M-band sur- 
veys towards young embedded stars all concluded that this as- 
signment is roughly consistent with the data. The ice mixtures 
found to provide the best fits in the literature dKerret alJll993l 
IChiaret alll994l Il995l Il998l iBoogert et all2002bl) are almost 
exclusively pure H2O-CO mixtures containing 25 or 5% CO, 
but with temperatures ranging from 10 to 100 K. Occasionally 
good fits are found with irradiated H2O-CH3OH mixtures. The 
identification of a phase of CO mixed in water ice is also indi- 
rectly supported by observations, since large columns of solid, 
generally amorphous, water ice are known to exist along the 
same lines of sight. Nonetheless, we find that the assignment of 
the red component to an amorphous CO-H2O mixture presents 
some inconsistencies between the observations and both old 
(high vacuum) and new (ultra-high vacuum) laboratory results. 

Laboratory spectra of CO co-deposited with water at low 
temperature (T < 20 K) show two distinct absorption peaks. 
One well defined peak near 2152 crrT 1 and one generally 
stronger peak near 2138 crrT 1 , although different experiments 
disagree somewhat on the peak position of the l ower frequency 
peak some placing it at 2136 - 2137 crrT 1 (San dfordet alJ 
198» ISchmitt et alJll989l) . This shift may be related to the de- 
gree to which the CO and the water are mixed, since CO de- 
posited on a water s urface shows absorpti on at 2138 cm" 1 . In 
recent work by (e.g.jDevlinll992tlPalumboll997[lManca et alJ 
2001 ; Collings et al. 2003) strong experimental evidence is pre- 
sented that the 2152 cnT 1 peak is due to the CO bonding with 
the hydrogen i n OH dangling groups as was suggested by e.g. 
Sch mitt et al.1 (119891) . In the same picture the 2137 crrT 1 peak 
is the normal CO bond site also found in pure CO but perturbed 
by a water ice surface. This holds in particular for CO trapped 
in micropores inside the amorphous water ice. Therefore the 
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2152 crrT 1 feature represents the principal binding site for CO 
on water, the presence of which would be strong evidence for 
CO interacting with an OH containing species. The 2152 cirT 1 
bond is suppressed by warm-up and is known to completely 
disappear for T > 80 K. Also irradiation by UV photons 
or bombardment by energetic particles destroys the bonds re- 
sponsible for the 2152 cuT 1 feature. In general, warm-up of 
a H2O - CO mixture also narrows the 2136 cirT 1 feature and 
shifts it towards the red. Irradiation also shifts the feature to 
the red, but in contrast to pure thermal processing broadens the 
band slightly. All of these processing effects are irreversible. 
Therefore the search for the 2152 cirT 1 feature is essential be- 
cause it has the potential to not only uniquely confirm the asso- 
ciation of the red component with OH-bearing species, but also 
to provide a sensitive temperature and processing indicator. 

None of our interstellar spectra show any evidence of the 
2152 cirT 1 OH-CO bond. The lower limits on the ratio between 
the optical depth of the red component and the 2152 crrT 1 fea- 
ture for the best quality spectra are shown in Fig.|^ They are 
compared to the temperature-d epend ent r atios from the labora- 
tory spectra bv lSandford et alJdl988l) and lSchmitt et al.l (lT989). 
Assuming that the laboratory spectra used are suitable as inter- 
stellar ice analogs, the limits on 8 sources are good enough 
to exclude ice temperatures below 60 K while an additional 6 
sources can exclude ice temperatures below 40 K. However, 
the profiles for ice mixtures at temperatures higher than 60 
K have FWHM of only ~ 7 cm 1 , making them inconsistent 
with the broader profiles observed along all lines of sight. It 
can be concluded that the available laboratory spectra in gen- 
eral exclude the low temperature non-processed hydrogenated 
ices, otherwise traditionally found to be consistent with inter- 
stellar spectra, due to the absence of CO in a OH bonding site 
(the 2152 cirT 1 feature). Alternatively, the interstellar ices may 
in general be strongly irradiated, but this requires an explana- 
tion why no other signatures of irradiation are seen along many 
lines of sight, such as OCN , aliphatics, etc. Strong irradiation 
may also be inconsistent with the strength of the typical UV- 
field in a dark cloud, since UV photons can only be formed as 
a secondary effect to cosmic ray hits due to extinction except at 
cloud surfaces. 

A significant problem in the interpretation of the red com- 
ponent is thus revealed. On the one hand, there is clearly evi- 
dence in the data for the presence of a CO component embed- 
ded in the water ice and for migration of a pure CO compo- 
nent into a porous water ice (Sec. 13. 6> . The shape of the red 
component, a high efficiency of CO migration and the general 
shape of water ice bands towards low-mass YSOs all suggest 
that most of the water ice column h as a fairly low tempera- 
ture jBooeert et ai]l2000bt [fh3|2002). On the other hand the 
absence of the 2152 crrT 1 feature seems to exclude this. It is 
a challenge for laboratory studies to explore how the CO-OH 
bonds can be avoided or efficiently destroyed under interstel- 
lar conditions. Possible explanations may include the destruc- 
tion of CO-OH bonds on interstellar time scales (> 10 3 years). 
Also, the formation of water molecules on grain surfaces may 
produce an ice structure different from structures obtained by 
deposition. Finally, the possibility that the red component is 
not due to CO interacting with water should not be ruled out, 
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Fig. 21. Lower limits on the ratio between the red component 
at 2136.5 cm" 1 and the 2152 cm" 1 OH-CO bond for the high- 
est quality interstellar spectra. The horizontal lines indicate 
the limits while the curves show the expected ratios from co- 
deposit ed mixtures of H?Q and CO. Solid line: H2O : CO = 
20 : 1 JSandfordet al.lll988h . Dashed line: H 2 : CO = 4 : 1 
dSchmittet all 19891) 

although it seems unlikely given the ubiquitous presence of wa- 
ter ice in grain mantles. 

For the red component, assuming grain shape effects are 
negligible, the column densities can be found from: 

N TC = 1 6 .0 cm" 1 x r max X A" 1 , (4) 

where A rc is the band strength of the bulk material. It assumed 
that the CO concentration small so grain shape effects can be 
ignored. The numerical factor is then the equivalent width of a 
t = 1 band. In the following a band strength for a water-rich 
mixtu re of A = 1.1 x 10~ 17 cm molec 1 is used JGerakines et alJ 
1199.4 . 

4.6. The blue component 

Since the blue component was only identified recently in the 
high resolution Keck -NIRSPEC spectrum of L1489 IR by 
Boo gert et alJll2002bl) . the only available statistics are from our 
work. In the present sample, it is found that the component is 
particularly prominent as a distinct shoulder in the profile in 
L1489 IR and in Reipurth 50. The feature is detected as a well- 
defined shoulder in many of the other sources such as HH 100, 
Elias 32, SVS 4-9 and IRS 63, but has a smaller ratio with the 

midd le componen t. 

lBoogertetal] (Eb02b) suggest that the blue component can 
be identified with CO mixed in a C0 2 -rich ice (C0 2 /CO > 1) 
or a mixture with less C0 2 but with significant amounts of 2 
and N 2 present. While reasonable fits can be obtained with a 
suitable C0 2 -containing laboratory mixture, we will argue in 
the following sections that another candidate explanation ex- 
ists, namely the longitudinal optical (LO) component from pure 
crystalline CO, which appears when the background infrared 
source is polarised. The presence of this component can be in- 
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Table 5. Polarisation fractions assuming LO-TO splitting 
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Fig. 22. The adopted model of the complex dielectric function 
of a-CO for p-polarised light. 



dependently tested by measuring the linear polarisation frac- 
tion at 4.7 fim. 

4.7. LO-TO splitting in crystalline a-CO? 

Laboratory spectroscopy of multi-layered crystalline CO (a- 
CO) using a p-polarised (light polarised parallel to the plane 
of incidence) infrared sour ce shows a splitting of the CO- 
stretc hing vibration mode (Chan g et alJ Il988t ICollings et alJ 
120031) . The splitting is due to vibrations perpendicular and 
parallel to the ice surface which are generally referred to 
as the longitudinal optical (LO) and the transversal optical 
(TO) modes, respectively. In this article model parameters for 
LO-TO spli t CO are derived using the laboratory work of 
IChangetalJ fl98R). unless otherwise stated. For thin ice layers 
(< 30 mono- layers), the two mod es show very narrow profiles. 
According to lChang et alJ (119881 the LO and TO modes have 
widths of 0.25 cm -1 and 0.85 cm -1 , respectively. The LO mode 
disappears for a single monolayer. For thicker layers (> 100 
monolayers), the TO mode exhibits a profile very similar to 
the profile seen in unpolarised light, i.e. cjq = 2138.5 crrT 1 
and FWHM = 2y = 1.5 cm -1 . The LO mode gives for 
thick layers a strong very narrow and blue-shifted profile with 
o)o = 2143 cm" 1 and 2y _ .5 cm -1 . Here €q = 1.84 is adopted 
as given in IZumofenl ll978). Unpolarised light will produce a 
single peak at 2138.5 cm" 1 very similar to that seen for amor- 
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phous pure CO ice. It is assumed that the dielectric functions 
for a-CO and amorphous CO are identical when subjected to 
unpolarised light. 

The dielectric function of a-CO for p-polar ised light is 
therefore modeled by two Lorentz oscillators ( iChang et alJ 
1988). The plasma frequencies of the two modes are assumed 
to be identical. A value of co p = 120 cm" 1 is found to repro- 
duce the strength o f both modes in the laboratory spectrum of 
IChang et alj fl988). This value is allowed to vary slightly to ac- 
count for uncertainties in the determination of the sample thick- 
ness. The LO-TO splitting model is corrected for grain-shape 
effects using the same CDE grains as for the middle compo- 
nent. After correcting for grain shape, io p = 140 cm" 1 gives a 
significantly better fit to all the astronomical spectra. The differ- 
ence compared to the laboratory data is considerable, also since 
a crystalline ice in space is not expected to have a different os- 
cillator density than the laboratory sample. Further laboratory 
experiments are needed to confirm if a discrepancy indeed is 
present. 

The necessary linear polarisation fraction at 4.7 /mi to ex- 
plain the blue component with the LO mode is found through: 



C 



total 



^unpolarised . -^-LSClt. p + (1 -flCuJ, 

1 - -T4.7 



(5) 



where P47 is the polarisation fraction at 4.7 fim while 
Cunpoiarised, ^LT.p and Clt.s are the absorption cross sections for 
unpolarised, p-polarised and s-polarised light, respectively, ft 
is a factor accounting for the fraction of the light seen as be- 
ing "p-polarised" by the CO molecules on a grain surface. In 
general, (3 — 1 /2 for randomly oriented particles, which can be 
realised by volume-integrating the projection on the plane of 
the incoming polarised light of the surface normal vector. If the 
grains are aligned along magnetic field lines, /3 can deviate con- 
siderably from 1 /2, but will never exceed unity. The CO ice is 
assumed to be entirely crystalline and the absorbing grains are 
assumed to be randomly oriented. It is important not to confuse 
the absorbing grains with the background grains producing the 
scattered and polarised light, since their properties may be dif- 
ferent. 
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Fig. 23. Spectra of sources with a clearly defined blue component. The Lorentz oscillator model for the middle component and 
the phenomenological red component have been subtracted. The plots on the left show the spectra compared to the grain-shape 
corrected laboratory spectrum of N2:02:C02:CO= 1 :5:0.5: 1 at 12 K by EAS. The plots on the right show the grain shape corrected 
LO-TO splitting model as described in the main text. All grain shape corrections use CDE grains. Positions of the gas-phase 12 CO 
v=l-0 (solid lines) and v=2-l (dashed lines) ro-vibrational transitions are also indicated, when present. 



The model fits are shown in Fig. [23] where they are 
compared to t he CO? -rich laboratory spectrum proposed by 
iBoogert et al.l (2002b) and EAS to fit the blue component, in 
particular in the case of L 1489. The derived polarisation frac- 
tions are given in Table[5] For L 1489 and Reipurth 50 the CO2 
laboratory profile indeed gives a good fit to the blue wing of the 
residual. For the other sources the fit is worse since the residual 
feature seems to be a factor of 2-3 narrower than the laboratory 
profile, although this may be an effect of line emission from hot 
CO gas as indicated in Fig. [53] In no case does the laboratory 
mixture give a good fit to the red side. It may be possible to 
construct a more complex ice mixture, which fits better. Also a 
better fit to the red side can be obtained by subtracting a mid- 
dle component of slightly smaller optical depth. However, since 
the slope of the blue wing is similar in all the sources, the exact 
same complicated mixture is required along all lines of sight, 
which seems very unlikely. 

On the other hand, the extremely narrow LO mode profile 
accurately reproduces the blue wing in the residual for all the 
sources. Taking into account the line emission from CO gas 
which clearly fills in the red side of the band in all sources 
but L 1489 and Reipurth 50, it is seen that also the LO-TO 
splitting scenario has problems reproducing the red wing. In 
L 1489 and Reipurth 50, excess absorption is clearly seen at 
2141 and 2142 crrT 1 , as previously mentioned. Nevertheless, 
the polarisation scenario has several advantages over the CO2- 
rich mixture scenario. First, since the polarisation fraction at 
4.7 yum is predicted, an independent test is available. L 1489 
and Reipurth 50 need very high polarisation fractions (> 50%) 
to explain the blue component. However, these two sources also 



have some of the highest K-bund polarisation fractions known 
in YSOs, as is indicated in Table |5] It is expected that the lin- 
ear polarisation fraction decreases with increasing wavelength 
if more of the light received originates directly from the emit- 
ting source rather than being scattered. In these sources, the 
polarisation fraction is seen to decrease from the H to the K- 
band and the maximum polarisati on fraction occurs o ffset from 
the mid-infrared position (Casali 1991; Wh itney et alJl997l) . It 
is not known if this behavior in general continues to 4.7 fim. 
If the grains are small and the emitting source is obscured at 
4.7 fim e.g. by a circumstellar disk, the polarisation fraction 
can be largely independent of wavelength up to 5 fim. Although 
few observations of polarisation fractions have been made at 
longer wavelengths than 2.2 //m, some protoplanetary nebulae 
have been found to show this behaviour, such as the Egg neb- 
ula, which has a constant linear polarisation fraction of 50% 
from 2.2 to 4.5 jum, thus re quiring grain sizes of less than 
~ 1 yum dKastner et aljl2002T) . Al so, the polarisation fractio n 
at 4.7 /mi for W 33A is known JChrvsostomou et aD fl996). 
The measured value of 10-13% is consistent with an upper 
limit of 15% from the decomposition of the ISO-SWS spec- 
trum. Curves of constant polarisation fraction have been drawn 
in Fig. El showing that the majority of the sources in the sam- 
ple will require polarisation fractions of less than 30%. 

The main problem in the determination of the polarisation 
fraction from the CO profile are the systematic uncertainties 
introduced by the assumptions on the fits of the middle and red 
components as well as the assumptions regarding the dielectric 
function of LO-TO split a-CO, especially concerning the width 
and strength of the LO component. 
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All the sources are thus consistent with a LO-TO split- 
ting scenario, but the result must be confirmed by imaging 
polarimetry at 4.7 /zm as well as further laboratory experi- 
ments using polarised light and high resolution spectrometers. 
It will have important consequences for the understanding of 
the processes governing the freeze-out of CO in the circumstel- 
lar medium if the blue component, and indeed also the middle 
component, is carried by crystalline CO. Note that the forma- 
tion of crystalline CO can occur at low temperatures, depend- 
ing on deposition rate and the adsorbin g surface and is not n ec- 
essarily a sign of processing of the ice ( Collings et al .120031) . 

5. Discussion 

5.1. Implications for the evolution and processing of 
CO-rich ices 

The most intriguing lesson learned from the presented data set 
is that the fundamental structure of the CO ice on grain man- 
tles seems to be practically identical along all observed lines of 
sight. It is well-established that the absorption profile of solid 
CO is sensitive to the conditions under which it is formed. This 
is an observation which suggests that the CO ices are formed 
under similar conditions regardless of the macroscopic context 
of the circumstellar medium or that the CO ice profile is invari- 
ant to differences which may exist in composition and grain 
shape distributions. 

It is well known that CO can efficiently freeze out on a 
given grain when it passes through a region of a dense cloud 
with the temperature and density most suited for the forma- 
tion of solid CO. In addition, an efficient method of desorp- 
tion at low temperature must exist, since CO is found in the 
gas phase in abund ances even in the dense st and coldest clouds 
(e.g. Ber gin et alJ l2001: Lad a et alJll994l) . An important des- 
orption mechanism seems to be via the release of chemical 
energy through reactio ns with rad icals produced by a cosmic 
ray induced UV field ( S hen et al although other mech- 

anisms connected with cosmic ray particles may contribute. At 
the same time ambient UV photons may only play a role on the 
surface of the cloud due to extinction. A given grain may also 
experience that all or parts of the ice mantle are evaporated and 
refrozen several times on the evolutionary pathway from dense 
cloud to protostellar envelope to circumstellar disk. 

In the light of this active processing which is expected to 
take place for any grain mantle, it is highly surprising that we 
see almost no change in the fundamental three profiles of solid 
CO for different lines of sight. The observational constraints 
presented here are thus that any processing, be it chemical or 
physical, must leave the profiles of the three components in- 
variant but may change the relative intensities. Typical lines of 
sight have 60-90% of the solid CO in a pure or nearly pure 
form. Models must be able to explain why this component of 
the CO has not experienced any kind of mixing with other 
species neither during formation nor as a result of subsequent 
processing. The results presented here thus seem to be consis- 
tent with a layered structure of the ices, which at least separate 
the CO ice from other species. Furthermore, the association of 
the red component with a simple FbO-rich ice mixture is prob- 



lematic. If the red component is indeed associated with H2O, 
it must be explained why it has the same profile for the quies- 
cent medium (e.g. CK 2), low mass YSOs (e.g. Elias 32), low 
mass YSOs in high mass star forming regions (e.g. TPSC 78), 
circumstellar disks (CRBR 2422.8) and high mass stars with 
abundant methanol (e.g. W 33A). Similar to the case for the 
middle component, this is unexpected since the profile of the 
red component should vary with the different temperatures and 
the different abundances of secondary species known to be as- 
sociated with the water ice, such as methanol and ammonia, 
in the varity of circumstellar and interstellar regions surveyed. 
These differences in temperatures, ice structure and composi- 
tion are observed through clear variations of the 3.08 fim H2O 
band and should be reflected in the shape of the red component 
as well. 

The implication is that solid CO seems to take part in little 
observable chemical processing around young stellar objects. 
Only the two sources L 1489 IRS and Reipurth 50 show some 
evidence for additional environments to the CO molecules 
through the excess absorption at 2141-2142 cirT 1 and even here 
the environments is not likely to have a very high column den- 
sity compared to the pure CO ice, since the features are so shal- 
low (t ~ 0.2). Chemical processing may still take place, but the 
products of any chemical reactions involving CO must be effi- 
ciently removed from the environment of the remaining CO. 

At the same time the variations in relative intensities of the 
three different components are consistent with physical pro- 
cessing affecting the total CO band profile, such as the evap- 
oration of the volatile pure middle CO component prior to a 
hydrogen-bonding red component as a result of thermal pro- 
cessing. This is supported by the observation that the red com- 
ponent is detected along every line of sight with solid CO. The 
presence of an exclusion region in Fig.Qolmay also support the 
mixing of pure CO into the water ice component as discussed 
in Sec.l3~6l 

5.2. Is the grain shape constrained? 

Even though the astronomical spectra can be so well fitted with 
the simple physical model presented here, it is a genuine worry 
that this is simply due to a degeneracy similar to that which 
haunts the mix-and-match approach. It is not possible to ex- 
clude a different combination of dielectric function and grain- 
shape model. It is, however, far easier to theoretically explain 
the common presence of pure CO rather than some specific and 
complicated mixture. Perhaps the most significant evidence for 
the applicability of a model using irregular grains simulated by 
a CDE shape distribution, is the evidence given when apply- 
ing the same combination of pure CO and grain shape correc- 
tion to o ther features. In particular it is shown bv lBoogert et alJ 
(2002b) and in this work that the stretching vibration profile 
13 CO along the line of sight towards two YSOs (one high mass 
and one low mass) is also consistent with pure CO indepen- 
dently of the grain shape. Furthermore, we have shown that the 
blue component, which is possibly the LO component of crys- 
talline CO, has a shape which is consistent with CDE grains as 
well. 
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Thus, under the assumption that the optical constants of the 
middle component are known, a strong constraint on the grain 
shape is given. The used CDE grains is not a realistic model for 
interstellar grains. The CDE model assumes grains of solid ice 
and includes grain shapes like thin needles and disks, which are 
not likely to be form ed. However, they are kno wn to simulate 
irregular grains well ( Bohren & Huffmanll 1983b and thus pro- 
vide a simple mathematical expression for this purpose. One 
possible scenario is that the pure CO exists as small irregu- 
lar clumps on top or inside the water ice mantles. This may 
unify the observational constraints given by the CO band with 
the constraints on the icy grain size distr ibution given by the 
3.08 fim water band JPartois et alJl2002l) . However, detailed 
calculations of the grain shape effects of simulated interstellar 
grains using for instance fractal models jFogel & Leungl 1997b 
are necessary to confirm the conjecture that CDE grains are 
appropriate for an irregularly shaped CO ice mantle. To fur- 
ther refine the adopted grain shape correction, additional solid 
state features from other species must be included. Ideally, it 
must be demanded that all observed solid state features are con- 
sistent with the same grain shape correction. However, since 
solid CO supplies one of the simplest and most easily observ- 
able bands, it likely provides the best single constraint. Other 
features which can constrain the grain shape correction fur- 
ther include the CO2 stretching vibration band at 4.27 fim and 
the CO2 bending mode at 15.2 fim. Also, having shown that a 
simple formula for the absorption cross section can reproduce 
the observed middle components, a convenient template is pro- 
vided against which more detailed models of grain shapes can 
be tested. 

5.3. Strategies for comparison with solid CO 
laboratory data 

The presented data provide new challenges for laboratory stud- 
ies of CO-rich ice mantles and some technical requirements 
become evident. The spectral features distinguishing the solid 
CO profiles of different sources show structures with widths 
less than 1 cirT 1 and much of the modeling requires high qual- 
ity laboratory spectra to investigate e.g. the structure of LO- 
TO split crystalline CO, which can have widths as small as 
0.2 crrT 1 . Accurate determinations of optical constants require 
laboratory spectra which are fully resolved. Thus, to match the 
present resolution of astronomical spectra, it is essential for the 
further study of interstellar solid CO to have laboratory data 
with a spectral resolving power of at least 0.1 cm -1 . 

The accurate astronomical CO ice profiles encourage a 
change of strategy when comparing observed solid state pro- 
files to laboratory data. As mentioned in Sec. 13.11 a mix-and- 
match approach is often employed to analyse the observed 
CO ice profiles. However, due to the many parameters govern- 
ing a given laboratory simulation (abundances of mixture con- 
stituents, deposition and annealing temperatures, irradiation 
parameters and the details of the laboratory setup) and given 
the low spectral resolution of the available laboratory spec- 
tra, degeneracies are introduced when trying to constrain the 
compositions of CO-rich ices in space using a mix-and-match 



approach. As has been shown in (e.g. iBoogert et al1 l2002b; 
Darto is et al . 2002) and this work, both the quality and quan- 
tity of astronomical near to mid-infrared spectroscopy is now 
sufficient that more advanced physical models can be applied 
to the shape of the observed solid state features. In particular, 
this requires on the laboratory side high resolution optical con- 
stants, but also detailed experimental and theoretical studies of 
the microscopic structure of simple ices under interstellar con- 
ditions. On the astrophysical side, it is evident that solid state 
features in general also probe non-chemical phenomena in ad- 
dition to the chemistry such as grain shape, size distributions 
and the overall temperature and density structure of the inter- 
stellar, circumstellar and disk material, which contain dust. The 
macroscopic properties of the objects observed are thus insepa- 
rable from the chemical properties of the ice mantles and must 
be modeled concurrently. Ideally, a single dust model includ- 
ing grain shapes and sizes, physical and chemical structure of 
the grain mantles and a macroscopic structure of the material 
along a given line of sight should be used to ensure consistency 
between the interpretations of different solid state features. 

6. Conclusions and future work 

Medium resolution M-band spectroscopy of a large sample of 
low mass young stellar objects has been used to explore the 
line profiles of the 4.67 fim stretching vibration mode of solid 
CO. A simple phenomenological three-component decomposi- 
tion reveals a remarkable similarity of the structure of the CO 
stretching vibration mode at 2140 cirT 1 between all observed 
lines of sight. A simple physical model is presented, which can 
reproduce the shape of all profiles to a high degree of accuracy. 

- The phenomenological decomposition of the profile is 
found to provide good fits for all observed sources. The 
decomposition uses three different components: a red 
Lorentzian profile centered on 2136.5 cm -1 , a middle 
Gaussian profile centered on 2139.9 cirT 1 and a blue 
Gaussian profile centered on 2143.7 cm -1 . Also the widths 
of the three components are kept constant such that only 
the relative intensities are varied between the sources. This 
approach reduces the information contained in a single line 
of sight to three linear parameters. 

- In 30 of the 39 sources a broad and shallow absorption band 
is detected centered in the spectral range 2165-2180 cm -1 . 
The central optical depth of the band is found to correlate 
with the optical depth of the red component of the CO ice 
band for low mass YSOs only. Higher mass sources tend 
to show excess absorption at 2165 crrT 1 . We propose that 
in addition to the traditional XCN band a second weak 
absorber centered at 2175 cirT 1 is present. Furthermore, 
due to the correlation with the red CO component, the 
2175 cirT 1 band may be due to CO in a new, unidentified 
binding site. 

- A simple physical model of the dielectric functions of the 
CO ice using Lorentz oscillators is described. It is shown 
that the middle component of the CO ice can be entirely ex- 
plained along all lines of sight by modeling pure CO with 
a single Lorentz oscillator and grain-shape correcting the 
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profile with a continuous distribution of ellipsoids (CDE). 
We suggest, that the CDE particles may work so well, be- 
cause they simulate irregular grain surfaces or small irregu- 
larly shaped CO clumps on top of the water ice mantle. The 
same model suggests that the different optical constants of 
pure CO obtained from laboratory experiments vary only 
due to a varying oscillator density (or porosity) of the ice. A 
variation of ice porosity in the lab can be caused by slightly 
different experimental setups. The observed middle compo- 
nents show no direct evidence for the presence of contami- 
nating species such as N2, O2 and CO2 mixed with the car- 
rier of the van-der-Waals interacting component and their 
concentration is likely less than 10%. 

- The blue component can be explained by the LO mode of 
pure a-CO, which appears when the background source is 
linearly polarised, although it cannot be ruled out that the 
component is carried by a C02-rich ice. If the blue compo- 
nent is due to the LO mode of crystalline CO, the polarisa- 
tion fraction at 4.7 fim is predicted in a number of sources. 
The derived polarisation fractions seem to correlate well 
with measured A'-band polarisation fractions from the lit- 
erature. The distribution of the relative contributions of the 
middle and blue components support the polarisation sce- 
nario. An important implication of LO-TO split CO is that 
a large fraction of the pure interstellar CO ice must be crys- 
talline. This would put strong constraints on the processes 
governing the formation and structure of solid CO in space. 

- The red component does not fit well with any simple 
H20-rich laboratory mixture due to the total absence of 
the 2152 cuT 1 feature, which is due to CO-OH bonds. 
Although thermal and energetic processing tend to destroy 
the CO-OH bonds, it also irreversibly changes the shape of 
the red component. Since CO adsorbed on many types of 
surfaces can produce a profile similar to the red component, 
it is suggested that alternative candidates which can explain 
the correlation with the 2175 cirT 1 are tested. Conversely, 
statistics of the relative contributions of the middle and red 
components support a scenario where pure CO migrates 
into a porous water ice upon warm-up. This scenario also 
predicts a strong 2152 crrr 1 feature. There is thus a signif- 
icant discrepancy between the data and all interpretations 
of the red component, which may require additional labo- 
ratory studies to solve. 

A number of observations and laboratory experiments are 
necessary to both confirm and to elaborate on the results pre- 
sented here. Further modeling will also significantly improve 
the understanding of solid state features of simple interstellar 
ices. 

- Modeling of the expected ice mantle structure using exist- 
ing dynamical models of collapsing protostars are neces- 
sary to theoretically understand the degree of mixing of the 
principal ice species which can be expected to exist at dif- 
ferent stages of the star formation process. Useful output 
from such a model would include the detailed composition 
of the ice as a function of mantle depth for a given grain. 
Clearly, such a model depends on the physical and chemical 
processes governing the adsorption of ice on surfaces, the 



surface chemistry of the grain, but also on the time the grain 
spends in different physical environments (density and tem- 
perature) during its lifetime. 

- Calculations of 'realistic' grain shape effects are required to 
explore which types of irregular grain surfaces can produce 
a grain shape correction similar to CDE particles. 

- UHV laboratory spectroscopy with resolving powers of at 
least 0.1 cirT 1 is needed to confirm the adopted dielectric 
functions for LO-TO split a-CO, to derive optical constants 
for other relevant CO-rich ices with narrow absorption pro- 
files and to further explore the effect of polarised light on 
astrophysical ices. 

- Future work should also include comparisons of other solid 
state bands of sources in the same sample such as H2O 
and CO2 with the results presented in this work. A physical 
model of the protostellar envelopes surrounding the sources 
using gas phase studies of rotational lines is required to 
fully understand the context of the observed ices and to sep- 
arate the possible contribution from foreground material. 
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Appendix A: Solid state line shapes 

We review here the physical argument for using Lorentz oscil- 
lators as a first order approximation to the dielectric functions 
of simple soli d state speci es. The formalism presented is mostly 
following lGadzukl dl987b . 

Common for both gas-phase and solid state lines is that any 
line shape from a spectroscopic transition is fundamentally the 
Fourier transform of the auto-correlation function of the rele- 
vant time dependent physical variable, q(t). 

1(a)) = I exp(-ia)(t-to))(q(to)q*(t-t ))dt. (A.l) 
*j 00 

This can be conveniently expressed with the Fourier trans- 
form Q(u>) = F{q(t), ai), via the convolution theorem by using 
that the auto-correlation function is the convolution of q(t) with 
itself: 

1(a)) = F(F~ l (F(q, to) x F*(q, a)))) 

= Q(a>) x Q(o)*. (A.2) 

The relevant solution for the pure harmonic oscillator is: 

q(t) = q exp(/w f - t]t), (A3) 

where a>o is the frequency of the oscillator and 77 is the dampen- 
ing factor determined by the finite decay time of the transition 
modeled by the oscillator. 
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It is well known that this solution, when inserted into eq. 
IA.2l vields a Lorentzian line profile, which collapses into a delta 
function when the oscillator is not damped, i.e. 77 = 0. Any 
model specific for solid state applications must take the inter- 
action of the oscillator with a background surface into account. 
A picture which is often used, is that the oscillator interacts 
via elastic collisions with a thermally fluctuating background 
of atoms or molecules. Every elastic collision will change the 
phase of the oscillator with a random phase shift. Such a phase 
shift destroys the correlation between the oscillator coordinate 
before and after the collision. If the collisions are random, but 
occur with a certain probability, an exponential correlation de- 
cay is introduced and the correlation function in ea. IA.2l is mod- 
ulated with a factor exp(-f/r c ), where r c is the average lifetime 
of the oscillator before a collision. It is straightforward to show 
that such a scenario will again result in a Lorentzian line shape 
broadened by l/r c . 

The effect is called de-phasing and is an often used concept 
in vibrational spectroscopy. Naturally higher order effects may 
produce more complex line shapes, of which the best known 
examples include non-homogenous broadening caused by non- 
random intermolecular distances. E.g. low surface coverage is 
known to have profound ef fects on the profile of CO adsorbed 
on metals dSomoriailll994l) . However, more advanced models 
may still be result in basic Lorentzian profi les mo dified with 
freque ncy-de pendent width functions as in iFanol J1961I) and 
iKubol H969). We conclude that the use of harmonic oscilla- 
tors and Lorentzian profiles provide a physically sound starting 
point for the theoretical study of the structure of interstellar ices 
using simple molecules such as CO as environmental probes. 
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Fig. B.l. 4 pm acquisition imaging of the SVS 4 cluster of 
young stellar objects in Serpens indicating the location of the 
two observed sources. The seeing is measured to 1.7 pixels = 
0"25. 



Appendix B: Notes on individual sources 

B.1. Ophiuchus 

IRS 43 is the closest source to the edge-on disk CRBR2422.8- 
3423 JThi et alJkOoll) . It has a reasonably deep ice band and 
shows broad ro-vibrational lines in emission from hot CO gas. 



lines blueshifted by ~ 100 km s 1 showing that the outflow ac- 
tivity is probably associated with this component. The depth of 
the CO bands are t = 2.1 and r = 1.2 for EC 90A and EC 90B 
north, respectively. This shows that the distribution of CO ice 
around young stars can vary significantly on physical scales of 
a few hundred AU. 



GSS 30 IRS1 is associated with a large well-studied reflection 
nebulosity. Almost no CO ice is detected along the line of sight 
towards this source although it is located in the most embed- 
ded regio n of the p Ophiuchi clouds according to the extinction 
map by Cam bresvl Jl999l) (A v > 10). The M-band spectrum 
of this source is distinguished from other sources by showing 
very strong ro-vibr ational lines in emission from gaseous CO 
dPontoppidan et al .120021) . 

B.2. Serpens 

EC 90 also known as CK 1 or SVS 20 is a 1 .6" binary young 
star, known to show signs of outflow activity l lHuard et alJ 
119971) . The CO ice band has previously been observed at low 
resolution bv lChiar et alJ i 19941) . but without resolving the bi- 
nary. Solid CO? (t = 2.6) at 4.27 yum has been detected by 
iGurtler et al.1 Jl996l) . We obtained well separated spectra of 
both components showing deep CO ice in both sources as well 
as cold CO gas. EC 90B additionally shows hot CO gas phase 



SVS 4-5 and 9: SVS 4 is a dense cluster of YSOs near the 
Serpens cloud core. The cluster is so dense, that it has often 
been confused with a single source. It was first resolved by 
lEiroa & Casalil dl989l) . who counted 1 1 bright sources in the 
cluster. Our acquisition image obtained under exceptional see- 
ing conditions (see Fig. IB. li confirms the number of sources 
and shows that no fainter sources are present to a limit of 
M4 n m = 11. Low resolution M-band spectra of the two br ight- 
est 2 yu m sources have been previously obtained bv lChiar et alJ 
dl994l) who observed SVS 4-9 and SVS 4-10, 5"to the north. 
We obtained simultaneous spectra of the two brightest 4 fim 
sources as indicated in Fig. IB. II the southern sourc e, SVS 4- 
9, being the same as observed by Chi ar et al 1 ill 994 while the 
northern source, SVS 4-5, is new. 

Deep CO ice is detected in both sources. The northern 
source has a typical narrow CO band, while the ice in the south- 
ern source shows a much broader feature consistent with the 
spectrum from the literature. 
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B.3. Chameleon 

These are the first ground-based observations of ices in the 
molecular cloud complex in Chameleon. The observations are 
difficult due to the extremely southern location of the cloud 
(6 ~ -11°). Consequently the sources never rise above an air- 
mass of 1 .7 as seen from Paranal, which makes proper correc- 
tion for telluric features difficult. Additionally, gas phase obser- 
vations are near impossible since the velocity shift compared to 
the telluric lines will be very small. 
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Cha INa 2 (source 1 in Table 2 of IPersi et all dl999l) ') is the 
only source with near-infrared excess detected by ISOCAM in 
the Chameleon I North a cloud. It has a class I SED and is 
surrounded by a small reflection nebula. It is estimated that the 
source is only extinct by 17 magnitudes (Pe rsi et alJl9 99). The 
source may be creating the large CO outflow observed in the 
area. This is one of the few sources which has a CO ice profile 
showing only a contribution from the red component. 

B.4. Corona Australis 

RCrA IRS 7 A and B There is some confusion in the literature 
regarding the identification of IRS 7 in the CrA cloud. At 2 and 
6 cm the cloud is domin ated by a stro ng binary radio source 
with a separation of 14'.' 2 jBrownl 119871) . It has been suggested 
that the two radio sources represents a strong bipolar outflow, 
since no obvious near-infrared counterparts to the radio sources 
were known. We detect two bright point sources at 4.8 fim lo- 
cated at the positions of the radio sources, and can therefore 
confirm that it is unlikely that the radio emission is from a 
bipolar outflo w from a single source. This was also found by 
IWilking et alJ 1 119971) who imaged the source at 10 /im and de- 
tected IRS 7A, but not IRS 7B, perhaps indicating that IRS 7B 
has a very deep 9.7 /im silicate absorptio n feat u re. 2 — 5 jjm 
spectra of IRS 7 was obtained by Tanaka et al. ( 1994), but it 
is unclear which source was observed due to the confusing K- 
band field. A narrow band M VLT-ISAAC image obtained at a 
seeing of 0735 is shown in Fig. IB. 21 

B.5. Orion 

TPSC 1 and TPSC 78 are two extremely red sources located 
in the trapezium cluster within 2'of the BN/KL region, with 
K-L > 6.1 and 4.85. respectivelv lEada et all ll200Ch . Their M- 
band spectra show very broad CO ice bands with the narrow 
component almost entirely missing as well as strong gas phase 
absorption. Also evident in both sources is strong absorption in 
the CN stretch region around 4.62 /mi. The ices are likely to 
have been affected by the close proximity to young high mass 
stars. 



Reipurth 50 IRS is an intermediate-mass source with a large 
associated reflection nebula similar to that of Cha IRN, GSS 30 
IRS 1 or L 1489. The s ource and the VLT-IS AAC spectroscopy 
is described in detail in Dartois et al. (2003). 




RCrA 



IRS7B 
I 



IBS 7 A 



Fig.B.2. Mnb image of the region around RCrA IRS 7 with the 
sources A and B indicated. Neither IRS 7A nor B is detected 
at 3.21 jjm due to the presence of a deep water band. IRS 7B 
is detected in the /T-band as a compact, but resolved reflec- 
tion nebula.The positions of the two sources are IRS 7A: a = 
19 01 55.16, 5 = -36 57 20.9 and IRS 7B: a = 19 01 56.25, 
6 = -36 57 26.9, J2000. 
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Table 3. Best fits to the CO ice band 



Cloud 


Source 


center(mc) 


T(rc) fl -'' 


r(mc)"- f 


T(bc) 1 ' 


Vel. corr. d 


CO gas vel.'' 


x 1 






cnT 1 








km s -1 


km s _1 




pOph 


IRS 42 


2139.78 +0.02 


0.072 ± 0.006 


0.53 ± 0.02 


0.05 ± 0.01 


+22 


+20 


0.82 




IRS 43 


2140.01 ±0.01 


0.15 ±0.02 


1.75 ±0.02 


0.07 ± 0.02 


-20 


_ 


1.76 




IRS 44 


2140.08 ± 0.08 


0.08 ± 0.03 


0.34 ± 0.05 


0.07 ± 0.03 


-20 


-12 


0.85 




IRS 46 


2140.24 ± 0.43 


0.09 ± 0.03 


0.11 ±0.06 


0.08 ± 0.06 


-20 


_ 


1.81 




IRS 48 


2139.47 ± 0.07 


0.01 ±0.01 


0.11 ±0.02 


0.03 ± 0.01 


-20 


_ 


2.36 




IRS 51 


2140.12 + 0.01 


0.42 ± 0.02 


4.54 ±0.15 


0.33 ± 0.03 


+23 


_ 


2.30 




IRS 54 


_ 


< 0.06 


< 0.12 


<0.06 


-20 


_ 


1.89 




IRS 63 


2140.01 +0.01 


0.28 ± 0.02 


1.69 ±0.03 


0.35 ± 0.02 


+23 


+ 18 


2.84 




WL6 


2139.94 + 0.03 


0.12 ±0.01 


2.04 ±0.12 


0.29 ± 0.05 


-20 


- 


1.45 




WL 12 


2139.82 + 0.03 


0.11 ±0.02 


0.79 ± 0.04 


0.07 ± 0.04 


-20 


-9 


2.10 




GSS30 IRS1 


2139.22 + 0.32 


0.06 ± 0.02 


0.07 ± 0.06 


<0.04 


+23 


+20 


0.57 




CRBR 2422.8 


2139.89 + 0.04 


0.47 ± 0.03 


4.38 ± 0.60'" 


0.13 ±0.07 


-17 


-9 


2.35 




Elias 32 


2139.95 +0.01 


0.39 ± 0.02 


1.76 ±0.05 


0.14 ±0.03 


-24 


— 


1.04 




VSSG 17 


2139.86 + 0.13 


0.28 ± 0.30 


1.04 ±0.36 


0.25 ± 0.32 


+22 


+20 


0.30 




RN0 91 


2140.25 +0.27 


0.43 ± 0.02 


0.25 ± 0.02 


0.15 ±0.04 


-19 


_ 


2.85 


Serpens 


EC90A 


2140.07 + 0.01 


0.24 ± 0.01 


2.09 ± 0.03 


0.06 ± 0.02 


-23 


-35 


9.84 




EC90B 


2140.00 + 0.01 


0.19 ±0.02 


1.11 ±0.02 


0.21 ± 0.03 


-23 


-35/-120 


5.12* 




EC 82 


2139.82 + 0.12 


< 0.03 


0.24 ± 0.05 


0.05 ± 0.05 


-22 


-18 


0.72 




EC 123 


2139.50 + 0.50 


< 0.12 


0.45 ±0.15 


< 0.09 


- 


- 


0.47 




CK2 


2140.06 + 0.05 


0.40 ± 0.07 


5.11 ±0.60 


0.23 ±0.12 


- 


- 


0.73 




SVS 4-9 


2139.93 +0.16 


1.07 ±0.10 


1.63 ±0.30 


0.42 ±0.14 


-22 


- 


0.52 




SVS 4-5 


2139.89 + 0.07 


0.59 ± 0.09 


3.32 ±0.71 


0.23 ±0.14 


-22 


- 


0.60 


Orion 


TPSC 78 


2140.13+0.04 


0.21 ±0.01 


0.12 ±0.01 


0.06 ± 0.01 


+8 


+8 


5.98 




TPSC 1 


2140.60 ± 0.88 


0.17 ±0.03 


0.06 ± 0.01 


0.09 ± 0.06 


+8 


+8 


0.74 




Reipurth 50 


2140.02 ±0.01 


0.412 ± 0.003 


0.51 ±0.01 


0.28 ± 0.003 


+7 


+ 10 


12.1 


CrA 


HH 100 IRS 


2140.01 ±0.01 


0.25 ±0.01 


1.32 ±0.01 


0.17 ±0.01 


-23 


— 


82.5* 




RCrA IRS7A 


2139.82 ±0.06 


0.39 ± 0.04 


1.41 ±0.17 


0.48 ±0.11 


-24 


-24 


1.14 




RCrA IRS7B 


2139.55 ± 0.07 


0.77 ± 0.06 


1.74 ±0.27 


0.56 ±0.10 


-24 


-24 


0.75 




RCrA IRS 5 A 


2140.08 ±0.02 


0.42 ± 0.02 


3.33 ±0.18 


0.53 ± 0.04 


-22 


- 


3.40 




RCrA IRS5B 


2139.98 ±0.10 


0.33 ± 0.06 


2.34 ± 0.57 


0.47 ±0.17 


-22 


- 


0.77 


Cha 


Cha INa 2 


— 


0.27 ± 0.02 


< 0.03 


< 0.06 


+3 





0.43 




ChaIRN 


2139.54 ±0.16 


0.17 ±0.05 


0.31 ±0.10 


0.05 ± 0.07 


+3 


- 


0.42 




Cha IRS 6A 


2139.82 ±0.05 


0.55 ± 0.06 


3.22 ± 0.47 


0.28 ±0.11 


+3 


- 


0.36 


Vela 


LLN 17 


2140.30 ±0.13 


0.32 ± 0.02 


0.19 ±0.03 


0.14 ±0.04 


-15 


- 


1.00 




LLN20 


2139.74 ±0.13 


0.29 ± 0.04 


0.38 ± 0.08 


0.12 ±0.07 


-10 


- 


2.05 




LLN 33 


2139.70 ±0.06 


0.33 ± 0.03 


0.94 ± 0.09 


0.21 ± 0.05 


-10 


- 


1.04 




LLN 39 


2139.97 ±0.10 


0.03 ±0.01 


0.02 ± 0.01 


0.004 ± 0.003 


-10 


- 


7.72 




LLN 47 


2138.78 ±0.43 


< 0.01 


< 0.02 


< 0.06 


-10 


-10 


3.24 


Taurus 


LDN 1489 IRS 


2140.04 ±0.04 


0.36 ± 0.02 


0.58 ± 0.03 


0.23 ± 0.04 


-11 


-10 


0.99 


Add. 


GL2136 / 


2139.50 ±0.06 


0.17 ±0.01 


0.15 ±0.02 


0.02 ± 0.02 






0.96 


sources 


NGC 7538 IRS 1' 


2140.04 ± 0.03 


0.10 ±0.01 


0.19 ±0.01 


0.10 ±0.01 






2.09 




RAFGL 7009S 7 


2140.19 ±0.04 


0.89 ± 0.02 


0.66 ± 0.06 


0.19 ±0.03 






3.82 




W33A' 


2140.90 ±0.08 


0.75 ± 0.05 


0.59 ±0.10 


<0.05 






1.01 




RAFGL 989^ 


2139.90 ±0.30 


0.18 ±0.03 


0.44 ±0.10 


0.14 ±0.03 






0.98 



" Solid CO column densities of the red component can be calculated using eq.|4| 
b Solid CO column densities of the middle component can be calculated using eq.|3| 
c Uncertainties are 3<r. 

d Calculated absolute velocity shift of the spectra, including Earth motion and LSR velocity of the parent molecular cloud. 

e Measured absolute velocity shift of CO ro- vibrational gas phase lines, where available. 

/ ISO-SWS archive data. 

" The optical depth given for the (saturated) CO band of CRBR 2422 is different from that o flThietalJi2003) since a laboratory profile 

was used rather than the emperical profile presented here. 
* The goodness-of-fit estimate for this source is affected by broad CO gas phase emission lines. 



